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iTils  report  discusses  the  design  end  synthesis  of  low-nass  and 
band-nass  one-dimensional  and  two-dimensional  soatial  filters  using 
arrays  cr  matrices  of  detectors.  The  synthesis  was  nerformed  for 
alternate  assumptions  of  negligible  end  significant  detector  gaps. 

Tne  detector  width  and  the  detector  weighting  functions  are  the  two 
optimized  parameters. 

A  freouency  analyzer,  using  a  narrow-band  one-dimensional  filter, 
was  desig:  and  simulated.  Tests  indicated  performance  as  expected. 
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I.  INTRODUCTION 
A.  GENERAL 

Recently,  the  detector  mosaic  concept  has  become  the  basis  of  many  new 
optical  processes.  In  long  wavelength  reconnaissance  systems,  terminal  guidance, 
and  pattern  detection  in  particular,  image  dissection  has  bee.,  found  a  pre¬ 
requisite.  These  systems  contain  a  re''atively  large  numb  r  of  input  detectors 
for  parallel  processing  and  recognition  of  the  object  with  h:‘gh  resolution  and 
high  value  of  signal  to  noise. 

A  spatial  filter,  as  defined  by  this  study,  also  uses  a  detector 
mosaic  as  the  input.  Its  function  in  space  is  analogous  to  classical  time 
filtering  in  that  it  scans  space  and  selectively  passes  spatial  frequencies. 

A  spatial  filter,  in  effect,  performs  a  convolution  in  space.  Its 
space  function,  or  transient  response,  analogous  to  the  system  fur.ction  in  time 
filtering,  is  constructed  by  weighting  each  detector  in  the  array-  The  weighted 
detector  signals  are  summed.  The  output  signal  of  the  summing  network  at  a 
fixed  position  of  the  detector  array  corresponds  to  one  point  of  the  convolution 
integral.  As  the  ar’^y  is  scanned  over  a  scene,  a  new  scene  is  generated  whose 
spatial  frequency  spectrum  is  the  product  of  the  filter  transfer  function  (the 
Fourier  transform  of  the  space  function)  and  the  original  scene's  spatial 
frequency  spectrum. 

Spatial  filters  can  be  one -dimensional,  in  which  case,  the  detectors 
are  rectangular  and  the  scan  motion  is  parallel  to  the  short  dimension.  Spatial 
filters  can  be  two-dimensional  in  which  case,  the  detector  dimensions  are  then 
determined  by  the  symmetry  of  the  filter  transfer  function.  For  example,  if  the 
transfer  function  is  circularly  symmetric  in  the  frequency  domain,  the  detectors 
also  will  be  circularly  symmetric  in  the  form  of  annuli.  A  two-dimensional 
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filter  can  take  any  form  so  long  as  it  has  qufldrantal  symmetry  with  respect  to 
both  frequencies'  axes;  that  is,  positive  and  negative  frequency  must  be  filtered 
identically. 

Ihe  spatial  filter  has  certain  fundamental  limitations.  For  example, 
a  high-pass  transfer  function  cannot  be  generated  because  it  would  require 
detectors  of  infinitesimal  size.  But  its  relative  simplicity  and  the  ease  with 
which  its  operation  can  be  generalized  (by  employing  adjustable  weights,  for 
example)  support  arguments  for  Investigation  of  the  spatial  filter  as  a  po¬ 
tentially  valuable  technique  in  electro-optical  processing. 

Previous  work  (Appendix  A)  has  reported  an  appix)ximation  method  for 
the  design  and  synthesis  of  a  one-diraensional  low-pass  and  a  band-pass  spatial 
filter.  Ihe  method  consisted  of  approximating  the  transient  response  of  the 
filter  with  a  linear,  finite-width  array  of  rectangular  detectors.  The  analysis 
ircluded  a  parametric  evaluation  of  approximation  accuracy  as  a  ftinction  of 
detector  size  and  number. 

The  continued  study  was  visualized  as  aiming  toward  a  synthesis 
procedure  which  could  be  useful  for  design  in  any  specific  problem.  To  demonstrate 
the  value  of  the  technique,  the  classical  problems  of  low  pass  and  band-pass 
filters  with  infinite  aotenuation  at  the  cutoff  frequencies  were  treated  as 
examples.  Ibe  eiYects  of  various  parameters  in  approximation  of  these  functions 
were  then  evaluated.  It  is  believed  that  practical  problems  do  not  present  such 
tight  specifications.  Therefore  the  approximation  in  a  real  case  would  be  per¬ 
formed  with  greater  accuracy. 

In  a  geneial  sense,  detector  size,  gaps  between  detectors,  and 
detector  geometry  are  parametrically  as  significant  in  spatial  filter  design, 
as  the  detector  wei^ts.  As  a  starting  point  in  the  general  analysis,  one  can 
assume  that  minimization  of  detector  number  is  important  and  that  close  approxi- 
macion  to  some  ideal  filter  response  is  desired.  Emphasis  on  various  aspects  of 
the  response  function  results  from  a  choice  of  approximation  criteria.  One  mi^t 
choose,  for  example,  to  optimize  for  least -mean -square -error  over  a  region  of 
interest,  or  to  optimize  for  least -mean -absolute -error.  Or,  as  in  time  filtering, 
maximally  flat  criteria  might  be  chosen.  To  a  limited  degree,  one  can  estimate 
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simply  the  errors  resulting  from  certain  approximation  methods.  Many  minimum-error 
optimizations,  however,  can  be  accomplished  only  throu^  iterative  computations 
best  implemented  with  a  hic^ -speed  digital  computer. 

Aside  from  its  value  in  spatial  filtering,  a  detector  array  has  the 
advantage  in  image  forming  of  mechanical  simplicity  over  a  single  scanned  detector. 
By  combining  this  stationary  image-forming  capability  with  the  spatial  filter 
concept,  recognizable  patterns  are  obtained  without  mechanical  scanning.  A 
frequency  analyzer  which  embodies  this  combination  was,  therefore,  proposed  and 
simulated. 

liius,  the  purpose  of  the  study  herein  reported  was  to  determine  in 
greater  derail  the  properties  and  parametric  potential  of  the  spatial  filter. 

Because  the  spatir.  filter  employs  a  detector  array  (either  one-  or  two-dimensional) 
or  variable  size,  and  because  the  fabrication  problems  of  detector  arrays  and 
processing  electronics  will  determine  the  feasibility  of  various  applications, 
establishing  the  parametric  sensitivity  of  th_  concept  is  a  matter  of  practical 
importance. 

The  initial  sections  of  the  report  develop  the  analytical  theory  for 
several  logical  approximations  in  both  one-  and  two -dimens ions.  Particular 
emphasis  is  placed  on  the  least -mean-square  error  (liBE)  approximation  because 
of  the  relative  ease  in  computing  the  optimization.  Several  analyses  develop 
the  theoretical  problems  of  special  geometries.  In  these  cases,  the  actual 
computation  of  filter  approximations  was  considered  unjustified  at  this  time 
and  the  analytical  results  are  summarized  in  order  that  the  study  of  these  cases 
may  be  continued  later  if  desired. 

The  later  sections  summarize  the  nximerical  and  experimental  work  and 
compare  the  various  filter  approximations.  The  practical  implications  of  the 
experimental  results  are  also  developed. 

B.  SUMMARY  OF  PROGRAM  GOALS 

Hie  program  started  on  1  July  1964  and  was  organized  into  the  following 

tasks : 

1.  Optimization  of  detector  size 

2.  Extension  to  two-dimensional  synthesis 
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5.  Background  discriinination  and  target  detection 

4.  Frequency  analyser 

5.  Multicolor  correlation  synthesis. 

A  mid-term  presentation  was  made  to  Mr.  P.  Nutting  of  ARPA  on  26 
October  196k  to  report  the  resultc  of  the  contract  and  to  concur  on  the  direction 
of  the  remainder  of  the  effort.  Mr.  Nutting  suggested  that  more  emphasis  be 
placed  in  the  detector  optimization,  considering  detector  gaps  as  another  parameter 
and  de-emphasize  the  color  correlation  and  target  detection.*  The  program  was 
scheduled  accordingly.  Figure  1  shows  the  schedule  and  performance  of  each  phase. 

1.  The  first  phase  entailed  the  evalxiation  of  an  optimization 
scheme  idiich  could  be  used  for  detector  width  and  vei^ting  function. 

2.  The  second  phase  was  an  extension  of  the  synthesis  procedure  to 
two-dimensional  detector  arrays.  A  least-mean-square -error  criteria  was  developed. 

3.  The  third  phase  consisted  of  a  s^Tithesis  procedure  for  practical 
detector  arrays  including  detector  gapis.  The  procedure  developed  in  the  f^  rst  and 
second  phase  was  implemented  in  this  phase. 

U.  The  fourth  phase  consisted  of  a  numerical  solution  of  the  above 
three  phases  to  evaluate  the  effect  of  detector  parameters  on  the  filter  response. 

5.  The  fifth  phase  consisted  of  developing  a  theory  design  and 
fabricating  a  frequency  analyzer  for  experimental  verification  of  the  spatial 
filter  theory. 


* 

The  importance  of  gaps  should  be  observed  in  practical  design.  In  a  one- 
dimensional  or  two-dimensional  mosaic  array,  fabrication  and  interconnection 
requirements  necessitate  some  gaps  between  detectors.  Therefore  a  practical 
synthesis  technique  should  recognize  their  existence.  In  addition,  con¬ 
siderations  of  economy  often  prescribe  a  maximum  number  of  detectors  for  a 
particular  application.  It  is  not  clear  that  a  uniform,  closely  packed 
arrangement  will  achieve  the  most  precise  synthesis.  Variable  sizes,  non- 
uniform  distributions  and  provision  of  gaps  may  be  preferable.  An  opti¬ 
mization  analysis  should  investigate  gap  width  as  a  parameter.  The 
numerical  solution  discussed  :n  this  report  shows  that  gaps  can  be 
introduced  to  advantage. 
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C.  SUMMARY  OF  PROGRAM  RESULTS 

In  sunimrizing  the  results  of  the  spatial  filter  study  ve  conclude: 

1.  Minimization  of  least-mean-square  error  is  a  successful  method 

of  approximating  si«tial  filter  fuiKtions.  For  many  approximations  the  conventional 
technique  coincides  with  this  approach. 

2.  Detector  gaps  with  regular  spacing  tend  to  degrade  the  filter 
perfonnance,  but  in  certain  casta  (Section  IV)  with  a  limited  number  of  detei^tors, 
variable  gaps  can  improve  the  overall  filter  response. 

Two-dimensional  square-response  filters  can  be  approximated 
with  a  cross  detector  array. 

When  a  square -response  filter  is  approximated  with  a  square 
detector  array ,  the  approximation  in  the  scop  band  imp'^ves  on  axes  between  the 
X  and  y  axes,  but  the  i»iss-band  ripple  is  worse. 

5-  Narrow -band  spatial  filters  can  be  satisfactorily  simulated 
with  transparencies  and  the  operation  of  a  non-scanning  frequency  analyzer  has 
been  denscnstratec . 

6,  One  additional  approximation  method,  detector  density  distri¬ 
bution  proportional  to  space  function  slope,  has  been  explored  analytically. 

7-  The  space  domain  least -mean-square-err:)r  (LMSE)  method  produces 
accurate  approximations  with  a  minimum  of  computation. 

These  results,  along  with  the  specific  parametric  data  are  significant 
in  any  future  filter  design. 
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SYN'ngSIS  OF  QNE»DIMENCIOML  SPATIAL  FILTER 
A.  I.EAST-MEAN -SQUARE -ERROR  APPROXIMATION 

1.  Lov-I^ss  Filter 

fbr  a  linear  array  of  detectors^  the  overall  detector  transfer 
f^in=tlon  Is  the  B«tien  of  the  Individual  detector  transfer  functions.  In  a 
unifora-slze-detector  array ,  the  transfer  function  of  the  1th  detector  having  a 

velght  and  centered  at  the  point,  j  1  -  |)  W  (where  W  is  the  detector  width) 
would  be  '  / 


(m)  .  A 


sinnnW 

nm 


exp 


2i-l 

2 


U  j  (2mn) 


(Il-l) 


'Ihe  transform  of  the  array  is  therefore 


+n 


2 


sinitmW 

jTin 


exp 


i=-n+l 


2i-l 

2 


WJ  (2.nm)j 


(II-2) 


Because  of  the  symmetry  of  the  transform  with  respect  to  positive  and  negative 
frequency,  Equation  (lI-2)  can  be  written  as 


F*  (m)  ==  2 

nm 


A. 

1 


cos 


[(2i-l 


(II-3) 


vUi  ideal  low-pass  filter  has  the  following  pass -band  and  stop- 
band  characteristics 


Fpj  (m)  w  1 

0  <  j  mj  <  j 

m 

c  1 

(II-4) 

^SE  “  ° 

|m|  >  |n._^ 

(II-5) 

where  m^  is  the  cutoff  frequency  of  the  filter. 


Page  II -1 


II  Synthesis  of  One -Dimensional  Spatial  Filter,  A  (cont. ) 


Report  No.  29B5 


Ihe  error,  iS,  in  the  pass-band  is 


^=F*~F 

0 

A 

m 

< 

m 

FB 

c 

{II-6) 


and  in  the  stop-band  it  is 


*'SB 

|m| 

> 

m 

c 

error,  E,  over  a  r^^gion. 

1 

A 

m  <  K, 

is 

(II -7) 


=  —  ( 
2K  J 


K 


dm 


-K 


Equation  (II-8)  takes  the  following  form 


(II-8) 


"  =  2r  <  f 


-m 


m 


n 


1  - 


sinmnW 

nm 


cos  {2i-l)nmW 


i=l 


dm 


.2  [ 

m 


sinnmW 

nm 


n 

z 

i=l 


cos  (2i-l)  ntuW 


dm)  (II-9) 


.0  have  the  Lf6E  we  require 


a  E 
a  A, 


=  0 


a  E 
a  A, 


=  0 


(Ti-10) 


a  E 
a  A 


=  0 
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These  conditions  produce  a  set  of  simultaneous  algebraic  equations  in  A^, 
^2  ^n  solved  for  these  parameters. 

Eqmtlon  (7T-9)  can  be  expanded  as  follows 


slnnmV 

nm 


n  -| 

^  A^  cos  (21“l)n!nW 
1=1 


slnnmW 

■2  •- 

nm 


^  ! 
y  cos  (2i-l)miiW  ^ 


K 

.2 

-  n  1 

1  r 

1 sinnmW 1 

)  A^  cos  (2i-l)nmW 

K  / 

1  ^  1 

m 

c 

-  i^l 

(Il-ll) 


As  the  set  of  partial  derivatives,  (lI-lO),  is  independent  of  the  order  of  inte¬ 
gration  in  (ll-ll),  the  set  of  Equation  (II-IO)  can  then  be  written  as 


^  £  1  r 

Ta^  -  K  J 


slnnmW 

nm 


2  ^  ^i  (2i-l)nmW  cos  (2i-l)nmW 


■2  (21-l)nmW  >  dm 


1 

^  K 


K 

r  / 

i  I 


sinnmW  ] 


n 

V 

-  i=l 


A^  cos  (2i-l)nmW  cos  (2i-l)nmW  dm  ( 11-12 ) 


After  simplification,  this  can  be  written 


Page  II -5 


A  (cunt.  ) 


Report  No.  2985 


II  Synthesis  of  One-Mmensionai  Spatial  Filter, 


1 

ic 


K 

.2 

n 

f  2 

sinnmW  1 

^  A^  cos  {2i-l)3nnW 

J 

nm  1 

0 

_  i-1 

cos 


(?i-l)ninW  dm 


2 

K 


sinnmW 

nm 


cos  (2i-l)nraW  dm 


0 


(II---' 5) 


The  set  of  Equations  (11-13 )  can  not  easily  be  solved  analytically.  A  numerical 
solution  is  clearly  indicated. 

If  the  limit  of  integration,  K,  is  allowed  to  become  very  large, 
the  numerical  manipulation  can  become  quite  complex.  For  practical  purposes, 
however,  it  is  possible  to  assign  an  upper  bound  to  K,  above  which  the  r'.,ponse 
is  not  of  any  interest.  This  mi^t  be  due,  for  example,  to  lack  of  information 
in  that  portion  of  the  frequency  spectrum. 

A  program  was  written  for  the  IM  7094  computer  for  the  cases 
of  n  =2,  4,  6,  8,  and  10,  and  K  =  4  m^.  The  results  of  the  numerical  solution 
are  discussed  in  Section  V. 

2.  Band-Pass  Filter 

Equation  (II-2)  applies  also  to  the  band-pass  case.  The 

coefficients.  A.,  however,  are  chosen  such  that  the  transform  satisfies  the 
1 

following  equations  instead  of  (II-4)  and  (II-5) 


F  (in)  =  Fpg  (  m)  =  1 


m 

c 


2 


m 


(  1 

F  (m)  =  Fgg  f  ra)  =  0  i 

r 

1 

< 

“c  ■  "2I 

i  r 

> 

m  /iffi 

c  +-2I 

(11-14) 


where  m^  is  the  center  frequency  of  the  band-pass  filter  and  is  the  band¬ 
width.  Applying  LMSE  criteria.  Equation  (lI-9)  -an  he  rewritten 
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E  = 


1  j 

r'"^  ■  =  1 

sinnmW 

2 

K  -  K  ,  > 

max  min 

c 

jtm  1 

n 


^  cos  {2i-l)nmW 
i=l 


dm 


1 


2an 

m  +  — 

c  2 


W 


i  ,  sinnmW 

I  ^  “ 


dm 
\  '  2 


nm 


A^  cos  (2i-l)nmW 


i=l 


dm 


r  max 

sinnmW 

2 

/  dm 

+  -X 

nm 

A^  cos  {2i-l)rtmW 


i=^ 


dm 


(11-15) 


Equation  {II-I5)  can  be  expanded  as  follows 


K 


E  = 


( 


max 


K  -  K  ,  j 
msx  min  ^ 


sinitmW 


Ttm 


K  . 
min 


n 


-I  2 


A^  cos  (2i-l)jtmW 


i=l 


dm 


(II-16) 


K  -  K  , 
max  min 


f 


n,  X  ^  f 

“c  ^  “  I 


<1-2 


sinnmW 


dm 

m - ^ 

c  2 


nm 


n 


)  A  cos  (2i-l)nmW 


L  i=l 


dm 


1 


For  LMSE,  the  set  of  Equations,  (ll-lu),  must  be  satisfied.  Differentiating 
+  :ie  integral,  we  can  write 


8  E  ^  1  / 

8  A,  K  -  K  .  j 
i  max  min  ^ 


f  max 


sinnmW 


nm 


rj} 

\ 


’min 


L  "1 

Li=l 


A.  cos  (2i-l)nraW 


cos  (2i-l)nmW  dm 


m  + 


dm 


re  2 

2 

j 

sinnmW 

K  -  K  , 

max  min 

f 

j  dm  ' 

m  -  -T- 

c  2 

1 

cos  (2i-i)nmW  dm 


(11-17) 


Equation  (ll-lj)  can  not  be  solved  analytically.  A  numerical  solution  can  be 


obtained  after  assigning  values  for  K  and  K  ,  . 

max  min 
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3 .  Lov-Fass  One-Difflensione  1  Spatial  Filter  Approxiniation 

IncludlnR  Detector  Gaps 

For  the  detector  in  the  original  one-dimensional  analysis  having 
width,  W,  we  substitute  a  new  detector  of  width,  <  W,  and  a  gap,  Wg,  (See 
Figure  ?.) 


Wg  -  g  g  <  1 

therefore  (ll-l8) 

W  = 
d  l^g 

W 

The  first  detector,  as  in  the  earlier  analysis,  is  still  centered  at  —  ,  the 

W 

second  at  5  g  > 

For  an  approximation  over  loops  and  using  2n  detectors 


—  =  nW  =  n  (l+g)  Wj 

2m  '  d 


(11-19) 


Equation  (II -3)  can  now  be  written 


F  (m)  =  2 


sinnmW 


d 


nm 


cos  (2i-l)nmW 


i-1 


(II-20) 


.  I  Win 
sin  I  i  m  T— '  ^ 

^  2  - - - cos  (2i-l)jtmW 

i=l 


nm 


for  loops  W  = 


2  m  n 


and 


siniim 


F  (ra)  =  2 


n 

n  m  (l+g) 


nm 


n 


i=l 


A^  cos 


(2i-l)  nm 


- 

2  n  m 

c 


(11-21) 
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Us.  Equation  (II-I3)  for  optimisation  on  the  basis  of  IKSE,  we  have 


sin  I  nm 


n  m  I  i-i-g 
c 


^  1 


A.  K 


''  A.  cos  (2i-l)jim  r— - 

1  2  n  m 


)s  [(2i-l)nni  —  dm  -  § 


2  n  m 


K  i 


.  ra  sin  nm  - - - 

r  c  n  (1+6. 


(II-22) 


cos  I  (2i-l)jim  I  dm  =  0 

L  ™c  J 


The  set  of  Equations  (II-22)  determines  the  coefficients.  A..  A  numerical 
solution  to  these  equations  has  been  made  on  the  7*^94  computer.  The  results  are 
discussed  in  Section  V. 

B  SLOPE  CRITERIA  FOR  DETECTOR  DENSITY 

In  addition  to  the  conventional  and  LMSE  methods  for  obtaining  detector 
veig.hts,  another  parameter,  detector  dv.nsity,  can  be  incorporated  in  the 
optimization  pixjcedure.  One  such  technique  would  be  to  have  the  detector  density- 
proportional  to  the  slope  of  the  ideal  space  function,  that  is,  the  detectors 
would  be  concentrated  \diere  the  function  changes  roost  rapidly. 


f(x)  = 


sir  2  n  m  X 
c 


(11-23) 


and  the  slope 


d  r(x, 


,  2  n  m  X  cos  2  n  m  x  -  sin  2  -  m  x 

1  c  c  c 


(11-24) 


The  detector  density  as  a  function  of  x  is  given  by 


N'(x)  =  u 


i  d  X 


(11-25) 
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vfhere  p  is  a  density  constant  determined  by  the  total  number  of  available  detectors 
and  the  truncation  point.  In  the  limit  as  the  truncation  point  and  the  number  of 
detectors  becomes  large. 


lim 

X,_  ^  CD 


f 


.00 


N’(x)  d  X  -  N, 


Total 


LiUl 

a  X 


d  X 


(11-26) 


Tlrerefore 


sin  2  n  m  x 


^Total  "  I  ■ 


X 


X. 


sin  2  n  ra  x 
c 


X 

sin  2  n  m  x  5 
c 


+  ... 


n  X 


X, 


(II-2T) 


sin  2  it  ra  x 


where  x,  ,  x, ,  x_,  etc.,  are  the  abscissas  of  the  - 

1^  3^  5^  "  n  X 


m'  'ima,  and 


sill  2  n  m  X 

Q 

X-,  X)  ,  Xr,  etc.,  are  the  abscissas  of  the  -  maxima. 

2^  ^  n  X 


Iherefore 


^Total  ^ 


t  f 


max  min 


(11-28) 


For  a  given  number  of  detectors, 


N 


p  = 


Total 


max  min 


(11-29) 


Ilie  increasing  series,  ^  |  f^^^  -  /  does  not  converge  because 


for  large  x^  it  becomes 


4  m 


n  2  n+1 


which  cleai'ly  diverges.  Hierefore,  in  a  practical  application,  some  finite  range 

for  X  must  be  chosen  and  /  { f  -  f  1  must  be  evaluated  in  that  region. 

!  max  minf  ^ 
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II- ‘  TWO-DIMENSIONAL  SPATIAL  FILTER  SYNTHRATR 

A,  CONV'ENTIONAL  APPROXIMATION  lECm’lQUE 

Assume  that  the  weights  for  a  two -dimensional  mosaic  of  detectors 
are  to  be  chosen  such  that  the  array  has  a  low-pass  characteristic  from  -m  to 

®c^  in  t  e  m^  direction,  and  from  -m  to  +m  in  the  direction,  as  shown 

-t-  ^  2  ^ 

in  Figure  3  i.e. , 


F  (m^,  ra^)  =  Fpg  (m^,  m^)  = 


F  (m^,  mp  =  Fgg  (m^,  m^)  .  0 


i  0  <  Imj^i  < 


rc 


0  < 


!  m 


^2' 


m-i  i  >•  I  m  j 

I  cp 

\  , 

^  i  m. 


V 


i! 

I  “‘2 


m 


'2 


(III-I) 


(III-2) 


The  synthesis  procedure  here  is  analogous  to  that  of  the  one-dimensional  filter. 
The  inverse  transform  of  this  ideal  filter  characteristic  can  be  written  as 


f(x,y)  = 


sin  2  n  m  x  sin  2  n  m  y 
1  P 


(in-5) 


f(x,y}  =  U  m  ffl 
c,  c 


sin  2  Ji  x  sin  2  n  m  y 

C^  c 

1 _ 

2  ji  m  X 


2  jt  m  y 
•^2 


(iii-iO 


The  two-dimensional  si»ce  function  f(x,y)  is  shown  in  Figx^re  3. 

For  a  mosaic  of  detectors,  the  approximation  can  be  carried  on  for 
2  loops,  3  loops,  etc.,  indepe.ndently  in  either  the  x  or  y  directions.  A  weight 
value  for  each  detector  can  be  obtained  by  either  the  conventional  method  of 

assigning  to  the  detector  the  -alue  of  f(x,y)  at  the  detector's  center  or  by 
applying  the  I/^SE  technique. 
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llie  value  of  f(x,y)  at  the  center  of  a  detector  located  in  the  ith  rov 
and  the  kth  column  is 


^ik  "  ^  ^ 


sin  2  !T  m  x.  sin  2  t  m  y, 
i  Cg  ''k 


='lS'k 


2  IT  X. 


2  n  y. 


(III-5) 


w. 


=  h 


sin  2  n  (2i-l)— nic^  sin  2  n  (2k-l) 


W 

2  r  (2i-l)  ~ 


Wo 

2  n  (2k-i)  ~ 


The  transform  of  a  i^ectangular  matrix  truncated  after  n  elements  in  the  x  direction 
and  p  elements  in  the  y  direction  would  therefore  be 


-  I  sin  IT  m,  W,  sin  n 

F  (m,,  m.)  -  4 - ^  - i-i 

1  2  2  !'  m,  W, 

IT  11 


"  “2  “2 


n 


I  I 

i=l  k-1 


sin(2i-l)jt  m  sin(2k-l)jT  W  m 

^  1  w  C 


(2i-l) 


(2k-l) 


(111-6) 

2  cos(2i-l)W^  IT  m^  cos(2k-l)W2  r 


m^ 


Equation  (ill -6)  can  also  be  written  as 


,  sin  n  m.,  w‘  sin  it  W  ^  f 

F  (m^,  m^)  —  n  m,  W,  n  m  V  /L  L  h^"^2i-l)n  (m^  +  m^) 

i=l  k=l  ^  ^ 


1  1 


sin(2i-i)rr  W  (m  -  m^  ) 
^  ^1  " 


sin(2k-l)n  (m^  +  m^)  +  sin{2k-l)iT  (m^  -  m^)! 


(III-T) 
1 


As  in  the  case  of  the  one -dimensional  filter,  F  m^)  can  be 

separated  into  two  factors 


0  (m^,  m^) 


sin  IT  sin  n  W, 

IT  m,  W,  IT  W„ 

1  1  2  2 


(111-8) 
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and 


n  p  sin 


f  (»i.  “2)  ^  ^ 


(2i-l)  n 


ffl 


cos 


(2i-l)  If 


i=l  k=l 


sin 


(2k -1)  W,  n 


la 


'2  J 


2i-l 


cos  (2k-l)  Tt 


(HI -9) 


ai^ 


2k -1 


ruction  and  n,)  are  called  the  envelope  function  and  edge 

function,  respectively. 


Scanning  the  two-dimensional  filter  in  the  x 


space  direction  corresponds 


response  function  on  the  axis;  i.e., 


to  filtering  with  the  one -dimensional 

=  0.  Ihe  envelope  faction  0  (^,  a^annlng  In  the  x  therefore 

sin  n  m  W 

takes  the  font  of  - H  1  ,.v,„ 

Tt  mf  Wf  '  the  one -dimensional  case. 

The  zeros  of  this  function  occur  when 


a. 


ra-.  W  =.*  a  or  at  m,  =  -i 

^7.  ■‘■1  1,^  W. 


(III-IO) 


where  is  an  integer. 

Scanning  the  detector  array  In  the  y  space  direction  corresponds  to  filtering 

vith  the  one-dimensional  response  function  on  the  axis,  I.e.,  m  .0.  For 

this  case,  the  envelope  function  reduces  to  ^ 


whose  zeros  occur  at 


where  is  an  integer. 


sin  T 
~^  “2  ^2 


m 


.  ^2 


(Ill-ll) 
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Scarjiing  at  some  arbitrary  angle,  9,  with  respect  to  the  x  axis  in 
tne  space  domain  corresponds  to  filtering  with  the  one-dimensional  response  function 
on  the  radial  line  making  an  angle  9  with  the  m^  axis.  For  this  line 


=  p  cos  9 
m^  =  P  sin  0 


(III-12) 


where  p  is  the  T-adial  distance  from  the  origin  in  the  frequency  plane.  Applying 
(lll-ll)  in  Equation  (lll-f)  we  find 


sin  (7;  p  cos  6)  sin  (n  W, 
^  p  cos  9 


p  sin  9) 


It  Wg  p  sin  9 


(in-15) 


Zeros  of  the  envelope  function  along  the  radial  axis  defined  by  9  are  therefore 

a. 


W,  cos  9 


=  P 


a. 


{lll-lk) 


sin  9 


=  7 


-1  '^l 

For  low  values  of  9,  i.e.,  9<  tan  ^  ,  the  zeros  of  sin  (n  W,  p  cos  9)  occur 

more  frequently  than  the  zeros  of  sin  (n  p  sin  9).  Figure  4a  shows  these  two 
terms  separately. 

As  the  value  of  9  increases ,  the  zeros  of  0^  wove  closer  to  the  origin 
and  the  zeros  of  0,  move  away  from  the  origin.  For  values  of  9  approaching 

-1  ^'1 

tan"  ”  the  negative  loops  of  0,  0^  become  increasingly  smaller  in  amplitude 
and  extent.  In  the  limit  when  p  cos  9  =  p  sin  9  the  negative  loops  in 
0  =  01  0^  disappear  completely.  A  more  quickly  attenuating  curve  as  shown  in 
Figure  4b  results,  where 


0  - 


si-,  (n  p  cos  9) 


^2 


n  p  cos  9 


(III-I5) 
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Returning  to  the  edge  function,  f  (m^ ,  ®p)>  on  the  axis,  •-  0, 
and  the  function  is  a  summation  of  cosine  terr.'»s  ^ose  zeros  are  (2i-l) 

(an;-!)  where  a  is  an  integer,  f  =  0)  has  a  maxima  at  -  ~  for  ^ 

even  values  of  a,  and  minima  at  ^  for  odd  values  of  Q.  liie  function,  F  (ra, , 

®2  =  0)>  has  therefore,  maxima  and  minima  on  either  side  of  each  —  which  is  a 

/  ^1 
zero  of  (m^,  ^  similar  result  exists  for  t.ie  m^  axis  and  . 

For  axes  at  an  arbitrary  9,  however,  the  edge  function  is  a  summation 
of  mixed  terms  such  as 

cos  f(21-l)  n  p  cos  9  j  cos  [{2k -l)  n  p  sin  9) 

For  small  values  of  9  where 

(21-1)  W  cos  9  >  >  (2k-l)  Gin  9 

/I 


the  zeros  of 

cos  j(2i-i)  •  p  cos  9j 

appear  more  often  than  the  zeros  of 

cos  [(2k-l)  TT  p  sin  91 
As  9  increases,  the  zeros  of 

cos  [(2i-l)  n  p  cos  9J 
move  to  the  left  and  those  of 

cos  [(2k-l  )  Wp  n  p  sin  9J 

move  to  the  right.  When 

(2i-l)  cos  9  =  (2k-l)  sin  9 


the  product  degenerates  to  a  cosine  squared  function. 
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B.  LMSE  APPROXIMATION  IK  FREQUENCY  DCMAIK 

In  the  liMSE  criteria,  the  weighting  function  of  a  detector  located  in 
the  i^^  and  column  would  be  a  value  A  ^  which  has  to  be  found  by  the  approxi- 
mating  procedure.  In  this  case,  formula  (III-6)  takes  the  form 


F  ( m-  ,  m  )  =  -^ 


n  p 

Sin  n  ra^  Sin  n  m^  ^  ^ 


m. 


m. 


Aik  ^'=*®[(2i-l)  i  ra^  W-, 


i=l  k=l 


(III -16) 


X  Cos  (2k-l)  ^  m^  I 


An  ideal  low-pass  filter  in  the  two-dimensional  case  has  the  characteristics 


PpB  mp  .  1 


0  <  m,  <  m 
1  c. 


< 

m 

1  2 1 

'  ^^2 

(III-17) 


and 


j  1 

"‘ll 

m  1 

"ll 

FgB  mj)  =  0  i 

\ 

1  ^2 

1 

i 

m  1 

"2' 

(III-18) 


where  m  and  m  are  the  cutoff  frequencies  of  the  filter  in  and  m^  directions, 

Cl  Cg 

respectively. 

The  error  in  the  pass -band  and  stop-band  is  similar  to  that  described  by  Equations 
(II-6)  and  (II-T).  The  total  error,  therefore,  in  this  case  would  be 

A 


E  = 


4  ^2 


A  d  d  m^ 


-'o  9o 


— ►  00 

Kg—#  OD 


(III-19) 
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£  = 


1 _ r^i  j  k  ^  ^  ^  ^2^2  ^  ^ 

^1^2  ill  li^  ^  ®p  ^  L 

“  i=l  k=i 


A.,  Cos 
ik 


I 

(2I-I)  IT  m..  W-  j  Cos 

X  X  i 


(2k-l)  n 


dm.,  dffl^  + 


.  j  I4,  ^  Sin  JT  Wg 

1=1  ic=l 


2  m, 

n  1 


A^k  Cos 


1  2  ^ 


{2I-I)  IT 


■"1  rh 


'm 

Cl  Cp 


Cos 


(2k-l)  IT  mg  Wg 


dm^  djUg 


(III-20) 


To  have  the  liCE  criteria  satisfied,  ve  require 


a  E 
Ta 


12 


=  0 


(III-2I) 


These  conditions  produce  a  set  of  np  simultaneous  algebraic  equations  in 
A^g,  ...  A^p  which  can  be  solved  for  these  parameters.  The  solution  of  these 
equations,  however,  can  not  be  easily  obtained  analytically.  E 'pending  these 
equations  similar  to  the  expansion  of  formula  (ll-ll)  through  (II-I3)  for  the 
one  dimensional  case  would  be  adaptable  to  a  compiter  i.nput.  A  numerical  solution 
for  cases  of  8  x  8,  10  x  10  and  20  x  20  square  matrix  is  obtained.  Ihe  results 
of  these  solutions  are  discussed  in  Section  V. 
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C.  LOW  BAND-PASS  OWO-DIMENSIOKAL  SPATIAL  FILTER  APPROXIMATION 
INCLUDING  DETECTOR  GAPS 


become^ 


For  the  two-dimensional  case  including  detector  gap,  Equation  (III-6) 


* 

F 


Sin  n  m,  W, 
1  d^ 

"i  "d. 


Sin  «  n 

2  S 


W. 

2  0 


z 

i=l 


)  A.,  Cos 

Lj  ik 


k=l 


“1  ‘'ll 


[11-20) 


:ii-2i) 


Cos 


2k-li  _  ,, 

— —  2  Ti  nig  Wg 


where 


Wi  -  (Ug^)  w 


and 


Wg  -  (l+gg) 


W  and  W,  are  the  width  of  the  detectors  in  the  x  and  y  directions, 
1  2 

respectively.  Detector  gaps  in  the  x  direction  and  y  direction  are  g,  W,  ,  and 

^  °1 

gp  W  ,  respectively.  An  optimum  approximation  to  the  filter  function  can  be 
2 

obtained  by  the  same  procedn.re  of  IJ^ISE  discussed  in  the  previous  sect  ' on.  The 
parameters,  and  must  be  chosen  for  the  approximation.  In  a  more  elaborate 
procedure,  one  may  optimize  not  only  on  the  basis  of  LMSE  for  the  weighting 
function,  but  also  for  and  g^. 


an 


[-7 
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IV,  SPATIAL  niEQUEtJCY  FILTER  DESIGN  BY  MINIMIZING  INTEGRAL  SQUARE 
A.  INTRODUCTION 

In  the  spatial  filter  design  methods  described  in  the  preceding  sections, 
detector  veightlng  coefficients  were  chosen  so  as  to  minimize  the  integral,  over  a 
prescribed  finite  region,  of  the  square  of  the  difference  between  the  desired  and 
approximate  frequency  domain  functions.  In  the  method  presented  here,  detector 
weigliting  coefficients  are  chosen  so  as  to  minimize  the  integral,  over  the  detector 
length  (area  in  the  two-dimensional  case),  of  the  square  of  the  difference  between 
the  desired  and  approximate  space  dcmain  functions.  But,  by  Parseval's  theorem  for 
Fourier  transforms,  the  latter  method  simultaneously  minimizes  the  integral,  over 
the  infinite  region,  of  the  square  of  the  difference  between  the  desired  and  approxi¬ 
mate  frequency  domain  functions. 

As  mentioned  in  Section  II,  the  minimization  of  error  in  the  spatial 
frequency  domain  was  over  a  limited  finite  region  because  extending  the  region  would 
have  led  to  computional  difficulties  on  the  IBM  709^.  The  method  presented  here 
overcomes  this  difficulty  in  that  it  minimizes  error  over  the  infinite  spatial 
frequency  region. 

Perhaps  the  main  advantage  of  the  method  of  detector  plane  fitting,  as 
shown  in  detail  below,  is  that  detector  weighting  coefficients,  for  all  of  the  cases 
treated  in  the  preceding  sections,  can  be  found  by  using  tabulated  functions,  thereby 
eliminating  the  digital  ccmputati ons . 

Furthermore,  as  shown  in  Sections  IV,D  and  IV, F,  detector  plane  fitting 
offers  a  simple  way  of  finding,  with  the  aid  of  a  digital  computer,  optimum  non- 
uniform  detector  widths  as  well  as  weighting  coefficients. 

To  see  how  Parseval's  theorem  applies  here,  let 

F{m)  =  F^^(m)  4  F^(m)  (iV-l) 

where  F(m)  represents  the  desired  spatial  frequency  characteristics,  F  (ra)  the 

a 

approximation,  and  F  (m)  the  error.  Also,  let  f(x),  f  (x)  and  f  (x)  be  the 

^  &  0 
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0 

corresponding  inverse  Fourier  transforms.  Then  let  the  integral  of  F  ‘"(m)  over 

6 

0  s  m  e»  defines  the  approximation  error  e  : 

«  =  f^F  ^(m)  dm  (lV-2) 

i  " 

Now  since  f(x)  is  real  and  thus  has  a  real  Fourier  transform,  l^rseval's  theorem 
can  be  written  as 


( *  /* 

I  f  ^  (x)  d  X  =  i  F  ^  (m)  dm  =  c  (lV-5) 

I  e  1  e 

■'o  0 

^  *  2 

This  equation  shows  that  if  a  design  procedure  minimizes  j  f^  (x)  d  x  by  fitting 

■6 

in  the  detector  plane,  i„  also  simultaneously  minimizes  the  spatial  frecuency 
approximation  error. 

In  the  above  equations  m  represents  the  spatial  frequency  variable  in 

cycles  per  unit  length.  But  in  the  remainder  of  this  section  U)  =  2m  in  radians 

per  unit  length  is  used  for  the  spatial  frequency  variable.  Although  this  change 

simplifies  notolion  here  by  eliminating  factors  of  2«  in  many  places,  it  requires 

chat  one  be  careful  in  relating  the  equations  of  this  section  to  those  appearing 

elsewhere  in  the  report.  Parseval's  theorem  can  be  written  in  this  way  because 

f  (x)  is  an  even  function  and  thus  has  a  real  Fourier  transform.  Frc^n  Equation 

(lV-2)  we  see  that  if  the  design  process  matches  f  (x)  to  f  (x)  by  minimizing 

& 


/  .f(x)  ~  f  vx)  ^  dx  (lV-4) 

"0  L  "" 

then  the  spatial  frequency  error  over  03m—®  is  simultaneously  minimized. 

Of  major  significance  here,  fitting  f  (x)  to  f  (x)  does  not  require 

& 

solving  a  set  of  simultaneous  equations  to  obtain  detector  weighting  coefficients. 
In  contrast,  the  present  ccmputer  optimization  solves  a  set  of  n  sim.ultaneous 
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equations  for  an  arraj  of  a  detectors;  and,  for  example,  in  tiie  case  of  n  -  tw¬ 
in  a  tv?o-diiriensional  ac’ray,  solving  the  equations  is  not  an  easy  task  even  lor 
the  IBM  709^. 

B.  QUALITY  OF  APPROXIMATION 


Because  the  ccmputer  optimization  minimizes  error  over  only  Osi  tt)i  nU5c 
■whereas  the  detector  plane  fitting  described  here  minimizes  erroi  over 
one  would  expect  a  smaller  error  over  O.'^uDirntuc  to  result  from  the  computer  opti¬ 
mization.  But  the  detector  plane  footing  described  in  Sections  IV,D  and  T/,?  leads 
to  detectors  of  variable  widt'" ;  ,  whereas  the  computer  optimizations  perfoimied  thus 
far  have  used  detectors  of  equal  widths.  It  will  not  be  too  surprising,  then,  to 
find  that,  as  a  result  of  the  latter  factor,  a  filter  designed  by  detector  plane 
fitting  may  have  less  error  than  by  frequency  dcanain  fitting.  Furthermore,  with 
equal  width  detectors,  the  fact  that  the  error  criteria  of  the  two  methods  beccme 
more  alike  as  n  a}„  increases  shows  that,  for  sufficient Ij-  large  n  resulting 

approximations  will  differ  negligibly. 


Although  they  will  not  be  used  elsewhere  in  this  report,  we  now  formulate 
expressions  of  approximation  errors  for  future  use  in  comparing  computer  oi't  ization 
■,ith  detector  plane  fitting.  Thus,  let  F  (tt))  represent  the  approximation  to  F  (u>) 
that  results  fron  frequency  donain  fitting  over  let  represent  the 

approximation  to  F({«)  that  results  ^rom  detector  plane  fitting.  For  an  ideal  low- 
pass  filter,  let  F(u))  =  1  over  ar.d  F(uj)  =  0  elsewhere. 


Then  F  ,(«))  minimizes 


ax 


■  ^  I  ^al 

0  ui^ 


(iv-5) 


and  F  ^(aj)  minimizes 

Utl 


= 


2  I 


T  2  '  w 

1  -  F  dti)  +  *  F  duj. 

^  j  -^5% 


(IV -6) 


Let  represent  the  errc.r  with  ¥^^(0)  over 


'0 


.2 

1  -  F  di'  + 

a2  ^ 


IVJ' 


c  ^ 


(rv-T) 
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er  i;esiRn, 


epuri*  riui. 


Ccmbining  integrands,  gives 


al-  ^aS> 


(lV-8) 


The  ratio  ^  thus  offers  a  measure  for  comparing  detector  plane  fitting  with 
frequency  plane  fitting  on  the  basis  of  integral  square  error  in  frequency  plune 
approximation . 

C.  OPTIMIMTION  OF  VEIGHTING  COEFFICIENTS  IN  A  OI^E-JlMENSIONAL  ARRAY 

OF  corfriGuous  detectors 


The  ideal  low-pass  spatial  frequency  function  again  serves  as  an  example. 
The  method  suggested  in  Section  IV, A  and  used  in  this  and  the  following  sections, 
applies  to  arbitrary  filter  functions;  the  details  should  be  easier  to  follow  for 
the  low-pass  case  because  of  its  simplicity  and  pr'’or  treatment  by  other  methods 
elseidiere  in  this  report . 


Thus,  consider  figure  6,  in  which  four  detectors,  each  of  length  , 

lie  in  the  Interval  0s:xS2n  (by  simmetry,  we  need  consider  only  the  half -plane 

X— 0)‘  As  suggested  in  Section  IV, A,  for  each  detector  the  integral,  taken  over 

the  length  of  the  detector,  of  the  square  of  the  difference  between  its  amplitude 
,  1  sin  X  ,  j  .  . 

^  X  minimized  with  respect  to  detector  amplitude- 


For  example,  in  the  second  interval 


'2  °  in/s'" 


-  Sg)  dx 


(lV-9) 


is  minimized  by  letting  a^  equal  the  value  found  in  solving  =  0.  Carrying  out 


the  minimization. 


t/2  ^ 


Then,  solving  =  0  gives 


■il. 


1  sin^x  ^^2  sin  x  2  , 

7  ”1 - r  -IT-  ^  ^2  ^ 

n  X 


sin  X 


=  ^  Si(n)  -  Si(|)l 


where  Si(x)  represents  the  s.ne  integral  of  x. 


(lV-10) 


(lV-11) 
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For  the  other  detectors^ 


^  jSl  (J)  -  0 


^3  =  ^  pi  (|-)  -  Si  {if)j  (l\^-12} 

%  ---  ^  |si  (2rt)  -  Si  {|i) 

D.  O^IMIZATION  of  VEIGHTING  coefficients  and  detector  WIDT  'S  for  a 
ONE-DB?EI'ISIONAL  ARRAY  OF  UNEQUALEY  SPACED  DETECTORS 

Consider  the  exp8.nded  sketch  in  Figure  7  of  the  interval  -  <  x  S  n 
where  the  diagonalli^  striped  areas  represent  the  integration  areas  of  the  error 
integral.  Assume  for  the  moment  that  a^  has  been  determined  as  above.  Now  consider 
what  happens ^to  the  error  integral  ;dien  the  second  detector,  which  originally  ex- 

tended  from  |  to  n,  is  replaced  by  one  that  extends  from  |  to  b  ,  thereby  forming 
a  gap  fran  b^  to  x. 

As  a  result  of  this  shortening,  the  entire  striped  aiea  no  longer 
represents  the  integration  area  of  the  error  integral;  to  the  right  of  x  =  b  ,  the 
dotted  area  replaces  the  striped  area.  Since  the  added  dotted  area  exceeds  the 
deleted  striped  area,  the  total  error  decreases. 

This  discussion  shows  roughly  how  shortening  a  detector  can  reduce 
total  error;  Equations  (IV-IT)  and  (r/-l8)  shovs  the  interdependence  between  b„ 
and  ag.  That  is,  the  procedure  in  the  previous  section  for  optimizing  a  no  longer 
applies  if  is  to  be  optimized  also. 

Further  consideration  of  this  shortening  process,  as  applied  to  the 
exasuple  of  four  detectors,  brings  out  some  additional  points.  First  of  all, 
shortening  the  detector  that  belongs  to  the  interval  0<x5|  will  not  reduce  the 
error.  Tnis^follows  fran  noticing  that  the  absolute  value  of  the  di inference 
between  — p-  and  a^  is  less  than  itself  throughout  the  interval.  Secondly, 

in  t,.e  intervral  x  :<■  x  S  ~  the  left  end  of  the  detector  gets  cut  off,  rather  than 
the  right  end . 
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The  procedure  given  here  applies  to  any  f(x)]  the  low-pass  filter 
prw.»ided  only  a  convenient  example.  Desigti  for  an  arbitrary  f(x)  can  start  by 
marking  approximating  rectangles  on  a  sketch  of  the  exact  function  prior  to 
making  tne  optimizing  calculations.  This  helps  to  organize  the  work  and  to  avoid 
imnecessary  calculations  in  regions  where  optima  do  not  exist. 

The  calculation  of  the  optimum  values  of  a^  and  in  the  above  example 
proceeds  by  first  writing  the  error  as 


^2  ^^2*  ^2^ 


=/;H 


\2 


r  It 


X  \  ;  ■  1 

1  ^ 


I  .  ^2 

sin  X 


X  i 


dx  (IV-I5) 


Then  ..'etting  ^  =  0  and  =  0  gives 


from  'idiich 


and 


rrcmi  which 


^*2 


i 


2  2  sin  X 


/2  " 


dx  +  2a, 


dx  =  0 


-'x/2 


2  51 


I 


2  sin  X  , 

/  -3^dx 


V,  " 

^2  ■  2 


siCb^)  -  Si(|) 
r. 

°2  ‘  2 


sin  b^ 


\2 


Ob,.,  n  b. 


-  a 


2; 

i 

I 


1 1  sin 
\ 


\  2 


=  0, 


(IV -14) 


(lV-15) 


(lV-16 


sin  b,, 
£  _ 2 

It 


(IV-IT) 
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Substituting  (iV-lj)  into  (IY-I5)  gives 


Si  (bj  -  Si  (i)  =  i2  -  ^ 


b. 


sir 


(lV-18) 


This  equation  is  readily  solved  vith  an  IBM  T09^  program.  For  rough 
calcui-atioriS  a  trial  and  error  procedure  is  uico  feasible;  substituting  the  value 
Oi  a^  that  uas  found  for  continguous  detectors  in  the  previous  section  into 


iln  b  a 


^  2 
o  ^ 

valued 


and  then  finding  b  in  tables  of  — furnishes  an  initial  trial 


E.  OPTIMIZATION  OF  VffilGHTLNG  COEFFICIENTS  IN  A  SQUARE  TWO-DIMENSIONAL 
DETECTOR  ARRAY 

The  two-dimensional  ideal  low-pass  filter  will  be  used  as  an  errample 
of  the  ootimization  method.  Let  ^  ‘  ^  ^ 


X 


0  elsewhere . 


Taking  the  Fourier  inverse  transform  gives 


(IV -19) 


for  the  exact  detector  function.  Assuming  that  the  detector  locations  and  areas 
have  been  prescribed,  the  error  for  each  detector  can  be  v-Titten  as 


'^ran  ”  j 


f 1 1  sin  X  sin  y  ^ 


nn 


dx  dy 


mn 

where  A^^  represents  the  area  of  the  mnth  detector. 

as 

Minlr,i2lr.g  by  setting  ^  .  o  gives 

mn 


(lV-20) 


'  2  sin  X  sin  v  ,  , 

^ ^  j  I - o  — Z —  — TT^’  dx  dy  +  2  a  A 

A 

mn 


(lV-21) 


P- ,  nr  ,  r\T  v 
‘■'•■Li'-  ■‘*“1 


X 
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Hence 


a 

mn 


X 


sin 


clx  dy 


(IV -22) 


Finally,  by  separating  the  double  integral. 


a 

mn 


Si  (x  )  -  Si  (x  ,)  Si(y  )  -  Si  (y  J  (iV-Fl) 

i  m'  rs-l'  ^•'n-1 


mn 


>diere  x  ,  x  .  , 
ra  m-i’ 


y^,  and  y^  ^  denote  the  boundaries  of  the  mnth  detector. 


F.  OPTIMIZATION  OF  ANNULAR  WIDTHS  AiH)  WEIGHTING  COEFFICIEtITS  IN  A  TWD- 
DBIENSIONAL  CIECULARL2f  SYMMETRIC  SPATIAL  IREQUENGY  FILTER 

The  analysis  in  this  section  uses  the  zeroth  order  Fourier -Bessel 
transform  (known  also  as  a  zeroth  order  Hankel  transform).  For  detector  weight¬ 
ing  lunctions  that  are  independent  of  polar  angle,  the  pertinent  transforms  becom,e 

..OB 

f(r)  =  I  ttF(uu)  J  (rj  )  duj  (lV-24) 

^0  ° 


and 

F('1})  =  ^  r  f  (r)  J  (r-ju)  dr  (lV-25) 

-0 

Here  we  apply  the  method  of  detector  plane  fitting  to  the  problem  of 
approximating  a  circularly  syrometric  ideal  low-pass  spatial  filter  by  a  set  of 
annular  detectors. 


Hence, 


F(«!)  =  1  for  uilltu 

c 


and 


F(-ju)  =  0 


for 

c 
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Taking  the  inverse  transform. 


-m  -5  1  w  u> 

f(r)  =  (u  J  (rtu)  doj  —  ts  {rm)  ~  J,  (r?o  )  (lV-26) 

j  o'  '  r  1  '  '  *  r  1'  c 

-^0  =-  ^  u 


Since  f(r)  appears  here  almost  in  the  form  of  the  function  that  describes 
Fraunhoffer  diffraction  at  a  circular  aperture,  we  may  take  advantage  of  this 
similarity  by  forming 

fir) 

,  2  ruj 
tt)  c 

c 

and  modifying  a  circular  aperture  diffraction  pattern  to  plot  vs  roj^  as  in 

Figure  8 .  ‘“c 

1.  Weighting  Coefficient  Optimization 


We  first  consider  the  optimization  of  weighting  coefficients  only, 
and  assume  that  the  detector  radii  have  already  been  chosen.  Figure  8  shows  the 
detector  configuration  and  the  desired  spatial  function. 

J^t  represent  the  error  over  the  area  of  the  innermost  detector, 
and  a^  the  v/eighting  coel’ficient .  Then 


=  2n 


r""!  f - 


[Ar)  -  aj 


r  dr 


(IV-2T) 


To  simplify  notation,  let  z  =  ruj^,  z^  =  ^  ~  — 2^  that  by  substituting 

(r/-26)  into  (r/-2T)  ^c 


2  I  1  I  ^2 

'l  “  lu^  \  —^  -a^'  z  dz 


Setting  =  0  gives 


(IV-28) 


2 


f  > 

1^-2  J  ^  dz  +  2aj^  I  ^  z  dz  I  = 


(lV-29) 


'M.  Born  and  E.  Wolf,  Principles  of  Optics,  Fergamon  Press  Ltd.,  1959, 
pp.  59^-596.  Page  IV -9 
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Hence, 


“ 


J^(  z  )d2 


Z 

1 


^  J  (z,  )  -  J  (0); 
L  o  1  o  j 

7..^  1 2 
1  ' 


(iv-50) 


Similarly,  for  a  detector  extending  frcci  r^  to 


a  = 
n 


J  (z  ,  )  -  j  (z  ) 
o  n-x'  o  n 

2  2/^ 

\  -  V-l 


(lV-51) 


where  a  =  — —  and  z  =  r  u; 

n  2  n  n  c 


2 .  Simultaneous  Optimization  of  Weighting  Coefficient  and  Radius 

The  procedxire  here  follows  that  shown  in  Section  IV, D  for  one- 
dimensional  detectors.  Using  the  preceding  example  of  four  annular  detectors,  ve 
shall  show  how  to  find  the  optimum  radius  along  with  the  optimum  weighting  coefficient 

Frcm  Figure  8,  ve  expect  that  reducing  the  outside  radius  of  the 

second  detector  should  reduce  its  error  in  approximating  '^l^^^c^  over  the  interval 

r'jtf 

c 

Letting  r  ^  denote  the  reduced  radius,  and  using  the  notation  on 
Equation  (IV-2P),  gives 


r^-2  Ij  (z)  -  f-2  jJ  {z)l 

-  ^1  -  0./  J  z  dz  (lV-52) 


^2  jj  (z)\ 


^€p 

As  in  Section  IV, D,  the  solution  follows  fran  setting  =  0  and  —  =  0 
Thus, 


^  =  2it  < -2  /  "^  J^(z)  dz  +  I  "z  dz  )  =  0 

'  I  ] 


(lV-53) 
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whence 


(r^-34) 


2  2  '  \ 

-  =  ^  ■  “2  ^ 


2 1  ^-2  ■  1  z  ,  ,  ^2  / 

\  -2  /  i 


(lV-35) 


’vAience 


2n  w  ll  \2 

-rf-  yj^.2^  -  “a  ^2  -  Ji  (^-2)  =  ° 


Q£  = 
2 


2  J,(z.2) 


(lV-36) 


(lV-37 


Equating  (iV^^k  and  (IV-5T)  gives 


1' 

2  2 
^-2  ■  ^1 


(lV-38) 


This  transcendental  equation  in  z_^  should  he  easy  to  solve  with  the  IBM  7094. 
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V. 


NUMERICAL  SOLUTION 


A .  INTRODUCTION 

To  evaluate  the  effect  of  finite  size  detectoro  and  the  effect  of 
detector  number,  or  truncation,  in  the  design  of  the  previously  described  low- 
pass,  band-pass,  one -dimensional  and  two-dimensional  filter,  numerical  solutions 
were  taken  from  an  IBM  709^  and  a  hybrid  computer  (HYDAC  2i+00),  To  normalisations 
vrei'e  performed  to  obtain  a  genei'alized  result.  First,  the  detector  width  was 
normalized  by  m^,  the  cutoff  frequency,  so  that 


m  W  = 
c 


2a 


(V-1) 


wiie^e  Q  IS  an  integer.  This  normalization  guarantees  an  integral  number  of  detector 
in  every  loop  of  the  envelope  function  . 

Second,  the  truncation  point  in  the  space  dcsnain  was  normalized  by  the 
cutoff  frequency  m^  i.e.. 


=  T  (V-2) 

where  is  the  abcissa  of  the  truncation  point  in  space  domain  and  t]  is  an  integer 
wi)ich  corresponds  to  the  number  of  loops  to  be  approximated  in  the  envelope  fimctlon 

The  IMSE  criteria  was  first  applied  to  an  array  without  detector  gaps 
in  order  to  ccsnpare  these  results  with  those  of  the  conventional  approach  as 
described  in  Appendix  A.  After  demonstrating  the  value  of  I^ISE,  this  method  alone 
W!is  used  for  approximation  of  linear  arrays  including  detector  gaps.  The  reason 
was  that  the  conventional  approach  could  not  ccanpensate  for  the  contribution  of  the 
gap  between  the  detectors. 

lor  practical  reasons,  LMSE  criteria  may  not  be  over  an  infinite 
trequcncy  range.  The  limit  frequency  band  chosen  for  the  evaluation  of  the  error 
depends  on  the  constraints  of  the  signal -and -noise  spectrum  in  the  stop-band.  In 
following  numerical  solutions,  the  approximation  has  been  optimized  up  to  four 
octaves  above  the  cutoff  frequency .  Approximation  for  a  different  frequency  band 
can  be  performed  in  the  same  manner. 
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B.  ONE~DIIffiNSIONAL  FIMER 

Figure  9  is  a  numerical  solution  for  a  low-pass  filler  trunCcated 
after  two  loops  of  the  envelope  function,  T1  =  and  employing  four  detectors. 
File  approximation  was  performed  by  both  the  conventional  method  and  the  IJ^BE 
metxiOd  with  ..ero  defector  gap.  As  shown  in  the  figure,  there  is  a  negligible 
eff^c^  on  both  the  p>ass-band  and  stop-band  response.  Furthermore,  the  weighting 

fun<..>ion  obtained  by  these  two  independent  methods  were  identical  out  to  three 
decimal  places . 


C.  TM3 -DIMENSIONAL  FILTER 

Numericai  solutions  for  a  two-dimensional  square  response  spatial 
filter  using  square  arrays  of  8  x  6,  10  x  10,  and  20  x  20  detectors  were  also 
obtained,  lo  generalize  the  solution,  the  frequency  spectrum  was  normalized  by 
®c^  and  as  in  the  case  of  one-dimensional  synthesis.  The  frequency  of  an 
8x8  aetector  matrix  along  the  x  axis  and  y  axis  Is  identical  with  the  one -dimensional 
low-p.iss  filter  (Figure  lO).  As  0,  the  scanning  angle  with  respect  to  the  X  axis, 
increases,  the  amplitude  of  the  harmonic  maxima  and  minima  in  the  stop -band  decrease 
as  illustrated  in  Figure  11.  The  same  phencmena  were  observed  for  the  cases  of 
10  X  10  and  20  x  20  detectors.  This  confirms  the  theoretical  discussion  of  Section 
III,B.  The  band  pass  ripple  increases  for  an  angle  0  between  0  and  90°,  reaches 
a  maximum  at  0  =  45  .  The  reason  for  this  increased  ripple  in  the  band -pass  may 
be  that  the  optimization  criteria  vhich  was  applied  for  0=0°  and  0  =  90°  only, 
thus  allowing  larger  ripple  for  0°'<,0*C90°. 

The  two-dimensional  approximation  was  perforaed  by: 

1.  Conventional  method  of  assigning  the  value  of  the  inverse 

transform  of  the  cube  at  the  center  of  each  detector  to  that 
detector;  and, 

2.  IMjE  to  minimize  error  out  to  4  m  . 

c 

The  weighting  functions  for  each  case  were  similar.  Differences  between 
the  weighting  functions  obtained  by  the  two  methods,  occurred  in  the  3rd  and  higher 
decimals.  The  weighting  functions  were,  of  course,  symmetrical  with  respect  to 
both  the  X  axis  and  y  axis.  Table  1  gives  the  value  of  these  functions. 
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A  Doint  of  interest  is  that  the  values  A,  .  for  high  i  and  j, 

i,J 

extremely  small.  The  value  of  A,  ,  for  8x8  detectors  is  about  O.O5  while  A,, 

yjP  XX 

is  0.81.  The  values  of  A,  ,  for  i  >3  and  j  ^5  are  all  less  than  0.01.  This 
indicates  hat  the  effect  of  desectors  located  beyond  the  3rd  row  and  3rh  column 
becomes  negligible.  Therefore,  one  can  use  a  reduced  raatrix  of  detectors  in  the 
form  of  a  cross  with  each  arm  of  the  cross  containing  four  rows  of  de-ectors. 


Figure  12  shows  a  tw^o -dimensional  spatial  filter  response  which  is 

sin  2«  m  x  sin  2n  m  y 
^1  ^2 

obtained  by  approximating  two  loops  of  the  - .  .  function 

^  Jtm  Jtm 

c, 

1 

with  a  20  X  20  detector  array.  A  comparison  between  this  curve  and  the  curve 
obtained  by  a  matrix  of  8  x  8  detectors  indicates  no  difference  in  the  band-pass 
response.  The  curve  indicates,  however,  a  better  approximation  in  the  stop-band  in 
the  region  m^  <“  m  than  in  the  case  of  8  x  8  detectors.  This  may  be  due  to 

the  motion  of  the  harmonic  maxima  and  minima  being  farther  from  the  cutoff  frequency. 

Attention  is  directed  to  the  large  ripple  in  the  stop-band  just  behond 
4  m  .  Had  the  approximation  been  carried  on  for  a  wider  band  than  4m  ,  this  ripple 
would  obviously  have  beer  attenuated  at  a  cost  of  an  increase  in  the  ripple  in 
either  band -pass  or  stop-band  before  4  m^. 

This  Iscussion  can  be  extrapolated  to  the  cases  of  circular  symmetry, 

low-pass  and  band -pass  filters.  A  series  of  circular  ring  detectors  centered  at 

Jl(r) 

the  origin  and  weighted  by  ———may  be  useu.  Since  the  first-order  Bessel's 

Jl(r) 

function  is  itself  damped  '/ith  increasing  r,  — —  decreases  more  rapidly  than 

.  Thus,  the  effect  of  detectors  located  far  from  the  origin  is  less  important, 
which  allows  earlier  truncation  in  the  approximation  procedure.  The  numerical 
solution  for  this  case  would  be  similar  to  that  of  the  one -dimensional  case  and, 
therefore,  was  nut  performed  at  this  time. 

D.  NUMERICAL  SOLUTION  INCLUDING  DETECTOR  GAP 

Figures  I5  and  l4  show  the  approximation  of  the  low-pass  filter  by 
using  an  array  of  r  lectors  with  varying  detector  gaps.  As  described  in  the 
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theoretical  discussion,  the  detectors  are  centered  at  the  same  point  as  in  the 
ca?.;  of  no  gap  The  vidth  ci  the  detector  is  decreased  v.o  allow  a  variahle  gap. 

The  phase  function  of  these  detectors,  therefore,  is  unchanged.  The  frequency 
specitrum  contributed  hy  each  detector,  however,  is  smaller  in  amplitude  and  more 
slosfly  damped  than  the  wider  detectors.  It  is,  therefore,  necessary  to  modify 
the  weighting  function  to  approximate  the  same  overall  spatial  filter  response. 

Due  to  this  general,  lower  amplitude,  the  weighting  function  of  each  detector 
must  he  increased  to  approximate  a  band -pass  characteristic.  Increasing  the 
weigiiting  function  also  increases  the  ripple  in  the  stop-hand  cesnpared  to  that 
obtained  with  the  wide  detectors.  These  results  were  confirmed  numerically.  Figure 
15  shows  the  effect  of  detector  gap  or  the  band-stop  ripple.  Frcm  this  figure, 
it  is  seen,  that  for  an  eight -detector  approximation,  the  band -pass  characteristic 
is  similar  to  that  with  no  detector  gap.  Ripple,  however,  is  developed  in  the  stop- 
band  region  of  U-5  m^,  8  m^  and  10  ra^.  If  less  ripple  is  required  in  the  stop-b-and, 
the  weighting  function  of  the  detector  .-.ust  be  reduced  accordingly.  Tnxi',  of  course, 
results  in  larger  ripple  in  the  pass -band.  The  slope  of  the  curve  aroui'.d  the  cutoff 
frequency  is  much  slower  because  o-^  the  narrow  detectors .  Figure  1^  shows  a  numerical 
solution  of  a  three -loop  approximation  for  various  detector  gaps.  The  fig’ire  illu¬ 
strates  closely  that  as  the  detector  gap  is  increased,  although  the  ripple  in  the 
stop-band  stays  fairly  constant,  the  approximation  in  the  pass  band  becemes  worse, 

A  cc»aparison  between  figures  indicates  a  sharper  slope  in  the  cutoff 
frequency  of  a  three-loop  approximation,  similar  to  the  case  of  no  detector  gap. 

The  case  with  detector  gaps  is  less  accurate  than  in  case  of  no  gaps.  A  better 
approximation  in  either  band-pass  or  stop-band  can  be  acccsnplished  by  weighting  the 
corresponding  error  in  the  LMSE  ccaaputation  more  heavily . 
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VI .  fhequeicy  analyzer 

A.  IIWRODUCTION 

I^ttern  ident-lfi cation  is  one  exaniple  of  a  promising  application  of  t-he 
=patial  filvcr  phenomena.  j.r  this  case,  a  criteria  for  categorizing  patterns  on 
the  jasic  of  their  spatial  characteristics  must  he  chosen.  One  of  the  simplest 
is  two-dimensional  power  spectra. 

There  are  a  number  of  methods  that  iw,ght  be  used  to  determine  two- 
dimensional  pwwer  spectra.  Such  spectra  are  obtained  conventionally  by  scanning 
a  slit  across  a  scene  and  analyzing  the  resulting  time- function  with  a  variable 

narrow-band  filter,  i.e.,  a  wave  analyzer;  or,  the  autocorrelation  function  is 
computed . 


An  alternate  procedure,  analogous  to  the  wave  analyzer,  would  require 
scanning  a  snccession  of  narrow-band  spatial  filters  across  the  scene  and  determine 
the  power  in  each  bandwidth.  Since  for  this  discussion  spatial  filters  are  weighted 
detector  arrays,  a  simpler  modification  would  use  a  single  array  and  variable 
weights.  Further  economic  argiiment  modifies  the  system  to  that  illurtrated  in 
Figure  I5.  In  this  case,  the  need  for  scanning  the  array  with  each  change  of 
weighting  function,  is  eliminated.  Instead,  the  equivalent  of  scanning  the  array 
is  accomplished  by  causing  the  weighting  function  to  travel  across  the  array. 

This  weighting  function  scan  can  be  achieved  practically  hy  using  a  tapped  delay 
line.  A  large  number  of  equivalent  scans,  therefore,  can  he  made  in  a  short  time. 


However,  as  will  be  shown,  the  methods  of  scanning  the  array  and 
scanning  the  weights  are  not  exactly  equivalent,  because  both  the  array  and  the 

weighting  function  are  truncated  or  framed.  Framing  introduces  a  distortion  into 
the  me' sured  spectra. 


B. 


DESCRIPTION  OF  SCANT®)  WEIGHT  FREQUENCY  ANALYZER 

Given  a  one -dimensional  weighting  function,  f  (x),  for  a  narrow-band 
spatial  filter,  having  a  center  frequency,  m  ,  and  bandwidth, 


2 

“  sin  (n^mx)  cos  (2n  m  x) 

.riX  0 


(VI-l) 
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We  initially  neglect  the  effect  of  a  finite  denslt;/  of  detectors j  that  is.  we 
aSuuiiiG  uhau  j.  (x)  is  approxiniated  perfectly  over  some  region  -  *  <x 
Whenever  j  x|>  the  truncated  f(x}  =  0 

j.he  effective  filter,  F  (m),  is  the  transform  of  the  truncated  f(x),  i.e. 


»  1  2  “ 
^  “  .J  (jt^x)  cos  (Pnm  x)  e 

3 


(VI-2) 


where 


-  F^(  m)  *  F  (  m) 


F^(  m)  s:  1  for 


G  2  c  2 


=  0  for 


m  I  <  m  - 


>  m  +  -— 
c  2 


sin  2n  m 

F  (  m)  =  - -_1 


In  examining  Equation  (VI-2),  we  note  that  the  convolution  of  F  (  m) 
ana  F^C  m)  degrades  the  edges  of  the  ideal  filter  function  and  produces  ripples  in 
hooh  t  -  pass-band  and  step-hand.  In  the  limit  as  becomes  very  large,  F  (  m) 
becomes  an  impulse  and  the  convolution  returns  the  ideal  filter  function. 

For  a  prac>.ical  situation  with  a  finite  detector  density,  F^(  m)  will 

net  oe  perxect,  xUt  will  exhibit  the  periodic  maxina  and  minima  demonstrated  in 
Appendix  A. 

In  add-ition  to  the  distortion  of  the  filter  characteristic  by  truncation, 
or  framing,  the  scene's  apparent  spectral  content  is  distorted  because  of  a  second 
xra.ue  imp^sea  oy  a  finite  focal  plane.  Assume  a  scene  having  a  -ne-dtmensional 
spectral  aistrihution,  G^(  m) .  If  the  scene  is  fra red  by  a  function,  G  (x),  where 
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Gg(x)  ^  1  for  ig  <  X  < 

=  0  for  ig  <  jx  I 

then  the  apparent  frequency  spectrum,  G*(  m) ,  the  the  convolution  of  G^(  m)  and 
Gg(  m) .  That  is 

G*(  m)  =  G^(  m)*  Gg(  m)  (VI-5) 

where 

m)  =  J-  e  to 

■‘■2 

=  ^  sin  S-mlj 


The  true  frequency  spectrum  of  the  scene  is  spread  by  G^f  m) .  Clearly 
Gg(  m)  sets  a  limit  on  both  the  resolution  of  spectral  measurement  and  on  the  lowest 
spatial  frequency  measurable.  In  other  words,  if  a  scene  contains  two  impulse 
functions  in  its  spectral  distribution  these  cannot  be  resolved  if  they  are  separated 
by  a  frequency  less  than  ,  because  the  frame  blurij'  the  two  impulse  functions  in 

the  apparent  spectral  distribution.  Furthermore,  frequencies  less  than  caruiot 

be  resolved  because  they  cannot  be  separated  from  the  impulse  always  present  in  dc 
(the  average  level  of  a  scene  must  be  positive). 

The  frequency  spectrum  of  the  signal  at  the  output  of  the  narrow-band 
filter  is  the  product  of  F*(  m),  the  effective  filter  function,  and  G*{  m),  the 
scene's  apparent  frequency  spectrum. 


the  weig 
fore,  is 
density, 
theri  is 


Notice  that  for  the  case  where  the  scene  is  scanned  across  the  array, 
nuing  function  is  framed,  but  the  scene  is  not.  Spectral  resolution,  there- 
limited  solely  by  the  truncated  weighting  function  and  finite  detector 
When  the  weights  are  scanned  and  the  scene  is  fixed  on  the  detectors 
the  additional  limitation  of  a  framed  scene. 
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In  the  following  simulation  experiment,  the  variable  transmissivity  of 
a  photographic  film  weighted  the  intensity  distribution  of  a  scene.  Mathematically, 
since  such  weighting  can  only  be  positive,  this  is  equivalent  to  adding  a  positive 
constant  to  the  weighting  function.  Comparing  with  Equation  (vi-2) 

F  -  (  ra)  =  r  r “  sin  (m^x)  cos  2nm  x  +  d]  ^ 

film  .  L  c  J 

"^1 


=  F*(  m)  + 


_D 

itm 


sin  2xm/^ 


(VI-U) 


The  constant,  D 

new  term,  ~  sin  2rtiat  is 
■  rtm  1 

frequency  characteristic. 


,  must  be  larger  than  2d,  m  to  keep  positive, 

centered  at  m  =  0,  and  it  distorts  the  filter's  low 


The 


C.  SIMULATION  EXPERIMENT 


The  experiment  simulating  the  frequency  analyzer  is  illustrated  by 
Figure  l6.  As  mentioned  above,  photographic  film  replaced  the  scanned  weighting 
function.  The  55-mm  film  was  driven  hy  a  slow-moving  transport  across  the  image 
plane  of  a  camera.  A  photomultiplier  behind  the  film  collected  the  light  from  the 
weighted  scene.  This  signal  was  plotted  by  a  pen-recorder. 

One  problem  in  the  simulation  was  construction  of  a  film  transmissivity 
proportional  to  the  desired  one -dimensional  narrow-band  characteristic.  Figure  IT 
Illustrates  the  procedure  used;  note  how  the  desired  function  is  plotted  in  the 
form  of  a  silhouette.  By  using  a  cylindrical  lens,  each  vertical  line  image 
appearing  at  the  film  plane  of  the  view  camera  has  an  intensity  proportional  to 
the  total  illuminatiun  of  the  corresponding  line  at  the  silhouette.  The  vertical 
cylindrical  lens,  in  other  words,  images  every  point  on  a  vertical  line  in  object 
space  to  every  point  in  a  vertical  line  in  image  space. 

A  negative  constructed  in  this  fashion  is  shown  in  Figure  l8.  This 
negative  was  then  used  to  construct  a  positive  transparency  on  jS-mm  film.  In 
processing  both  the  negative  and  positive  it  was  necessary  to  develop  to  a  7  =  1 
to  maintain  linearity. 


Page  VI -4 


VI  Frequency  Analyzer,  C  (cont.) 


Report  No.  2985 


The  figures  illustrating  the  film  processing  show  that  the  filter  chosen 
jxtended  only  to  the  first  minimuTn  of  the  envelope,  that  is,  to  where 


,  -  n  and  X,  =  ” 
1  1  zan 


Also  fan  was  chosen  such  that 


^  =  0.1 

m 

c 


D  2dim 


Referring  to  Equations  (VI-2)  and  VI-4) 


F_(  m)  =  ~  sin  20n 
2  Jtm  m 


(VI-8) 


,w.  ^  0  •  2ni 

F  ,  (  m)  =  F  (  m)  +  — - — -  sin  20n  — 

film  «m  m 


(VI-9) 


0.2m 


The  term. 


c  .  ni 

—  sin  20n  —  ,  is  negligibly  small.  Therefore 
'c 


(VI-10) 


Also,  since  Xp  -  (because  both  are  the  length  of  a  55”®m  frame). 


G„(  m)  =  — -  Sin  20n  — 
2  Jtra  m 

c 


(VI-11) 


m  m 

Three  "picket  fence"  scenes  with  frequencies  m  ,  ,  and  -§■  ,  and  the 

c  £-  5 

resultant  signals  when  scanoed  with  the  film  band-pass  filter  are  shown  in  Figures 
19>  20,  and  21.  The  effect  of  framing  the  scenes  J  Gq(  m)  is  seen  in  the 


Page  VI-' 


. . . . . . .  . . . . . . . . .  . .  . .  "  '"I . I"'  . . . . . . . . .  . 


VI  Frequency  Analyser,  C  (cont.) 


Report  No .  29-35 


gradioal  increase  and  decrease  of  the  trace  at  the  beginning  and  end  of  the  scan. 
Theoretically  a  picket  fence  of  frequency,  ,  contains  no  energy  in  its  second 

€m 

harmonic,  m^.  However,  the  cumulative  errors  in  the  spectral  measurement  produce 
a  small  signal  at  m  in  Figure  20. 

m 

A  picket  fence  of  frequency,  — ,  should  contain  a  signal  at  m  with  an 
amplitude  one  third  that  of  the  fundamental.  Figures  19  and  21  illustrate  the 
3:1  amplitude  ratio  of  the  m^  components. 

The  high  frequency  “hash"  at  the  edges  of  the  traces  are  presumed  due 
to  the  effect  of  transparent  tape  attaching  opaque  backing  at  the  edges  of  the 
filter  function. 


In  the  film  simulation,  the  following  suggestions  may  be  considered  for 

future  work: 

1.  As  noted  earlier,  by  exposing  the  film  with  too  much  contrast 
the  film  is  driven  beyond  the  linear  region.  A  reduction  in  contrast  would  make 
a  measurable  improvement  in  the  filter  function. 

2.  A  more  exact  simulation  could  be  made  by  imitating  the  finite 
detector  density  with  steps  in  the  silhouette  function. 

5.  Hy  splicing  together  various  strips  of  film,  it  should  be  possible 
to  enlarge  the  length  of  the  filter  function,  i.e.,  increase  ly  This  also  shotld 
make  it  possible  to  reduce  e:~^raneous  effects  such  as  the  very  high  frequency 
noise. 

4.  A  complete  set  of  bandpass  filters  and  test  patterns  might  be 
constructed  to  explore  completely  the  potential  of  the  frequency  analyzer  principle. 

It  is  also  possible  to  use  two  separate  films  in  a  beam  splitter,  with 
one  beam  for  positive  weights  and  one  for  negative  weights.  The  outputs  of  the  two 
are  subtracted  to  get  an  exact  equivalent  weighting  function,  thus  eliminating  the 
d-c  term  in  the  original  function. 
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VII.  ADAPTIVE  mmmum 

As  described  in  the  report,  there  are  numerous  parameters  which  affect  the 
performance  of  a  spatial  filter  including  detector  width,  number  of  detectors, 
detector  gaps,  detector  density  function,  and  optical  system  transfer  f* motion. 

As  the  analysis  of  these  parameters  combined  is  cumbersome,  this  study  was  directed 
toward  evaluawing  the  effect  of  each  parameter  independently  in  the  overall  filter 
response.  An  optimum  design  for  a  specific  filter  function  should  enjoy  the 
flexibility  of  varying  these  parameters  simultaneously. 

The  design  of  a  practical  detector  array  may  become  more  complex,  however, 
if  the  non-uniformities  of  detector  response  and  telescope  transfer  function  as 
well  as  detector  width  and  gap  tolerances  are  to  be  considered  functions  of  the 
f  cal  plane  dimensions.  These  difficulties  and  attendant  problems  can  be  avoided 
or  minimized  by  applying  the  adaptive  learning  technique  in  the  design  of  these 
filters.  This  scheme  would  entail  placing  an  array  of  detectors  in  the  focal  plane 
of  a  telescope.  Each  detector  signal  would  be  fed  to  a  variable  weighting  element 
which  would  be  adjusted  by  the  signal  of  a  comparator  network.  The  comparator 
would  sense  the  difference  between  the  actual  output  of  the  detector  array  and  a 
described  output  for  a  prescribed  input  function.  The  weighing  element  would  be 
adjusted  repeatedly  until  the  difference  becomes  zero.  At  this  time,  the  weighting 
function  values  would  be  fixed.  A  sequential  and  iterative  process  thus  could  be 
employed  to  realize  a  filter  with  prescribed  frequency  response  characteristics. 

Future  work  should  include  the  feasibility  evaluation  of  adaptive  learning 
techniques  for  realization  of  the  spatial  filters. 
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2rm  (RADIANS  PER  UNIT  LENGTH)  WHERE 
m=  SPATIAL  FREQUENCY  (CYCLES  PER 
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IDEAL  LOW-PASS  SPATIAL  FREQUENCY  FILTER  AND  A  DETECTOR  DISTRIBUTION 
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A  COMPARISON  BETWEEN  FREQUENCY  RESPONSES  OF  A  ONE- 
FILTER  OBTAINED  BY  CONVENTIONAL  AND  L.M.S.E.  METHOD 


THE  RESPONSE  OF  A  TWO-DIMENSIONAL  SPATIAL  FILTER  ON 


Report  No.  2985 


Figure  11 


w  ONJ^DJMENSON^  DETECTOR  ARRAYS  =2 
(FTITED  TO  RECTANGULAR  FREQ.  CHARACTERISTIC) 


Report  Nc.  2985 


Figure  :i2 


FREQUENCY,  CYCLES  /  UNIT  LENGTH 


Report  No.  2985 


Figure  I5 


TIMING  ^  1 

SIGNAL  S  SPECTRAL  POWER  FUNCTION  OUT 

INTEGRATOR  ° 


Figure  l6 


Report  No.  ?9S5 


Figure  17 


VIEW 

CAMERA 


-DIMENSIONAL  BAND-PASS  TRANSPARENCY 


Report  No.  2985 


Figure  18 
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APPENDIX  A 

SPATIAL  FILTER  SYNTHESIS* 


Bv  K,  Seyrafi  and.  G.  A.  Davison 


ABSTRACT 


^  ^slgn  and  synthesis  of  one-dimensional  low-pass  and 

snace^d^Rin^t^^  J  -Iters  are  developed  by  approximating  the 

^  response  of  the  filters  with  rectangular 

ors.  The  transforms  of  the  approximating  functions  are 
er  yed  an^ytically.  Numerical  solutions  to  the  function  are 

for^wafo^Phl^  ^  computer.  The  results  indicate  that 

i  tolerances  in  the  pass-band  or  stop-band,  a  small 
number  of  detectors  is  sufficient. 


presented  at  the  10th  National  IRIS.  Ihe  work  was  per- 
ormed  uiider  a  subcontract  from  General  Precision,  Inc.,  Libra- 
scope  Division,  and  supported  by  AFSC  Space  Systems  Division 
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I.  INTRODUCTION 

Once  it  is  given,  the  transfer  function  of  an  electrical  network  can  generally 
be  realized  by  a  proper  combination  of  lumped  passive  eles.ents  (resistors,  capaci¬ 
tors,  6Uid  inductors)  and  active  elements  such  as  vacuum  tubes  and  transistors  (Ref¬ 
erences  1,  2  and  5)*  The  exact  configuration  and  number  of  necessary  elements  is 
determined  by  the  required  pass-band  and  stop-band  characteristics  (ecg.,  ripple 
and  rate  of  attenuation).  These  requirements  have  led  to  systematic  design  pro¬ 
cedures,  such  as  Chebycheff  (equal -ripple),  and  Butterworth  (msximally  flat)  filters. 
There  are  sdso  systematic  procedures  for  realizing  filter  designs. 

A  similar  design  and  synthesis  has  not  yet  been  developed  for  space  filter¬ 
ing,  althou^  the  application  of  two-dimensional  Foiurier  transforms  to  spatial  fil¬ 
tering  has  been  studied  in  seme  detail,  and  there  has  been  empirical  design  in  low- 
and  band-pass  spatial  filters  (Reference  4),  Uiere  is  no  simple  analogy  between 
the  concept  of  electrical  filter  elements  and  synthesis,  and  space  filters,  such 
that  one  is  able  to  extend  the  well -systematized  techniques  of  electrical-circuit 
synthesis  to  spatial- filter  synthesis.  However,  one  can  apply  the  approximation 
techniques  usea  to  find  the  transfer  function  of  an  electrical  filter  \dien  the  im¬ 
pulse  response  is  known.  The  time  impulse  response  can  be  approximated  by  a  series 
of  trianguleir  or  rectangular  pulses  (References  5  and.  6).  The  transfer  function 
is  found  by  summing  the  transfer  functions  of  the  elemental  pulses,  each  with  its 
corresponding  delay. 

The  following  discussion  employs  a  suggestion  of  J,  A.  Jamieson  to  approxi¬ 
mate  the  transfer  function  of  a  spatial  filter  using  an  enseatole  of  rectangular  de¬ 
tectors,  and  it  develops  a  method  similar  to  the  tine-function  approximation.  A 
systematic  scheme  is  thereby  proposed  for  designing  and  synthesizing  low-pass  and 
band-pass  spatial  filters  using  a  variable  number  of  resolution  elements  or  letec- 
ors.  The  outputs  of  the  detectors  are  weighted  according  to  their  position  in  the 
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image  plane,  emd  are  then  summed.  It  is  shown  that  by  increasing  the  number  of 
detectors  or  resolution  elements,  one  cein  approximate  the  desired  transfer  function 
with  arbitrary  accuracy,  which  is  analogous  to  increasing  the  number  of  elements 
in  an  electrical  filter. 

XI.  PRINCIFLE  OF  APPROACH 

A.  LOW  PASS 

1.  Derivation 

We  begin  by  examining  an  ideal  one-dimensional  low-pass  spatial 
filter  with  a  flat  response  of  unity  out  to  the  frequency  m  =  cycles/unit  length 
and  thereafter  a  zero  response,  as  shown  in  Figure  lA.  The  inverse  Fourier  trans¬ 
form  of  this  function,  which  is  the  spatied  impulse  response,  has  the  characteris¬ 
tic  shown  in  Figure  IB. 

The  curve,  f(x),  of  Figtire  IB  can  be  approximated  by  an  infinite 
set  of  rectangular  pulses  with  a  uniform  width,  W.  A  summation  of  the  trsmsforms 
of  the  rectangiXLar  pulses  is  an  approximation  to  F(jm)  and  in  the  limit  as  W  — f  0, 
is  equal  to  it.  The  realization  of  each  pulse  is  simply  a  detector  of  width  W  and 
a  sensitivity  proportional  to  the  pulse  height  (Figure  2). 

Ihe  amplitude  of  the  ith  pulse,  A^  (Figure  2),  is: 

,  sin(2i-l)  ~  2  7tra 

A.  =  -  - -  for  i  an  integer  (l) 

1  n  W 


f^(x)  where 


Therefore,  the  pulse  can  be  represented  in  the  space  dcanain  by 

.  2i-l 
sin 

■d*)  =  I 


^  W2  nm 
2  c 


(2i-l)  I 


(i-l)  W<  X  <  iW 


(2) 


X  <  (i-l)  W 
X  >  iW 
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The  approximation  of  f(x}  will  therefore  be; 
i  =  00 


f(x)^ 


(5) 


i  =  -  00 


Practically,  it  is  not  possible  to  incorporate  an  infinite  number  of  detectors . 

But  because  of  nhe  attenuating  effect  of  the  ^  factor  for  hi^  values  of  x, 
the  ampuitwde  of  the  pxilses  beccanes  increasingly  smaller  and  smaller,  such  that 
elimination  of  those  detectors  has  sli^t  effect  on  the  trans:fer  function.  There¬ 
fore,  if  we  truncate  the  approx  'ation  after  n  deto  tors  ar^d  denote  the  truncation 
function  by  f*  (x),  we  have: 

i  =  +n 


f  (x)  =  )  f^(x) 


i  =  1-n 


The  Fourier  transform  of  f  (x),  namely  F  (jm)  will  be  the  sum  of  the  individual 
transforms  of  f^(x): 


+n 


F  (Jm)  = 


i=l-n 


sin 

2i-l  o 

2it 

L  2  c  J 

(21-1)  I 


.w 


sin  2nm  I  _j^(:a^)W 


2jtm  - 


(5) 


As  F»(jm)  is  an  even  faction  of  m,  the  exponential  factor  can  be  reduced  to  a 
cosine  function. 
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^  .  2i-l 

E  sin  — 

I  -TO 


Wffi  2it 

c 


W 


sin  2jnn  ^ 
nm 


i  =  1 


2i-l 

2  cos  W  27011 


(b) 


By  rearranging  the  terms  of  Equation  (6),  ve  can  write 

r2i-l 


n 


—  sin  nm  W 
It  TonW 


sin 


(^)W  2Ttm^ 


cos 


i=l 


“(21-1) 


(2|ll)  y  2itm 


(7) 


since 


sin 


f  W  2ian 

1 

COS 

w  2ion 

sin 

(2i-l)nW  (m^+m)J 

+  sin 

‘ 

(2i-l)itW  (ffl  -m) 
L  C 

[2  cj 

!  2 

Equation  (?)  can  also  be  written 


(8) 


as 


i  Sin  nmW  Y"  '2i-l)  rtW  (m^+m)  j  +  sin  f  (2i-l)  ttV 

It  innW  (ITTll - 


i=l 


itW  (m  -m) 
c  . 


(9) 


which  is  a  parametric  approximation  to  an  ideal  one -dimensional  low-pass  spatial 
filter. 

2.  Theoretical  Discussio- 

The  truncated,  finite-detector,  low-pass  equation  has  several 
interesting  properties  that  can  be  predicted  analytically  and  are  confirmed  by  tne 
numerical  results.  In  the  limit  as  tne  truncation  point  moves  out  to  infinity  and 

as  the  detector  size  approaches  zero,  the  space  filter  function  approaches  tne  ideal 
low-pass  case. 
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The  approximate  filter  function.  Equation  (j),  can  be  separated 

into  two  factors, 

sin  rtffiW 

^  ~  lonW  '  which  we  call  the  envelope  function,  and 


i  =  1 

edge  function. 

It  can  be  seen  that  the  envelope  function  (F  gure  ^A)  will  cause 
the  filter  function  to  go  to  zero  at  every  frequency  m  =  ^  for  which  a  is  an  integer 
The  edge  ^.unction,  4^(ffl),  which  is  a  summation  of  odd  cosines,  is  shown  in  Figure 
5B,  The  edge  function  has  a  maximum  at  m  =  |  for  even  vaJues  of  a,  and  a  minimum 
at  2  for  odd  values  of  a.  The  fultev  fur.ction,  F*(jm)  =  |  5Z!{in)  H^(m) ,  will  there¬ 
fore  have  a  maximum  and  minim-un  on  either  side  of  each  “These  harmonic  maxima 
and  miniiua  can  be  seen  later  in  the  numericel  curves.  The  effect  of  using  increas¬ 
ingly  smaller  detectors  (decreasing  W)  should  therefore  be  to  move  the  harmonic 
pairs  of  maxima  and  minima  farther  apart. 

As  n  is  increased  and  additioneQ.  odd  cosine  functions  are  in¬ 
cluded  in  the  edge  function  summation,  the  slope  of  the  filter  function  will  break 

more  rapidly  before  and  after  each  narmonic  pair.  Indeed,  it  is  easily  seen  that 
the  summation, 

n 

sin  (2i-l)  wjan 

c 

i=l 

is  very  similar  to  the  Fourier  eummatic»i  for  a  square  wave  having  a  period,  -, 

^  sin  (2i-l)  I  COB  (2i-l)  Wnm 
i=l 


cos(2i-l)  V^.-an 
2i-l 


sin  (2i-l)  cos  (2i-l)  Wnm 


2i»l 


,  which  we  call  the 
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For  those  situations  where  Wm  =0.5^  the  curves  are  periodic 
square  wave  approximations  modified  by  the^  —  factor  and  limited  by  n  terras. 

For  the  curves  where  Wm^  'C  O.5,  the  coefficients  are  no  longer  exact.  The  effect 
of  extending  the  truncation  point  should  be  therefore  an  increasingly  shai'per  slope 
at  the  cutoff  frequency,  ra^,  and  also  at  the  harmonic  maxima-rainima. 

From  Figure  IB  it  ceui  be  seen  that  the  maximum  usable  detector 
width,  W,  would  be  ^  larger  detector  width  could  not  approximate  the  periodic 

sin  2jtm  x 

nature  of  the - —  function. 

X 


B.  BAND  PASS 

1.  Derivation 

The  mathematics  developed  in  the  previous  section  was  for  a  low- 
pass  filter.  A  band-pass  filter  can  be  obtained  from  the  difference  of  two  low- 
pass  filters. 


For  a  band-pass  filter  with  a  center  frequency,  m  ,  and  band- 

c 

width,  we  can  assume  that  the  filter  is  the  subtraction  of  two  low-pass  filters  - 

A|n 

one  with  cutoff  frequency  of  m  - ^  and  the  other  with  cutoff  frequency  of 

+  “2*  principle  of  superposition  of  a  function  and  its  transform,  we  can 

write  the  inverse  transform  of  the  band-pass  filter  as: 


1 

.in  2it(m^  +  ~)  X 

sin  2n(m 

c 

-  2  ;  X 

_  ^ 

2  cos  2nm^x  sin  nfmx 

X 

X 

nx 

- 

- 

(10) 


f(x) 


2  sin  roCanx 
jtx 


COS  2jan  X  =  - 

c  n 


sin  Ttfimx 
-inix 


COS  2jta  X 

c 


(1.1) 


2.  Theoretical  Discussion 

Equation  (11)  reduces  to  a  low-pass  filter  for  a  center  frequency, 
=  0,  and  will  have^"J^^^  characteristics.  When  —  is  large  compared  to  1,  f(x) 
behaves  like  a  cosine  function  (cos  2m^x)  damped  by  a  modul  '■  Ing  —  ^  ^  function. 
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As  the  bandwidth  of  the  filter  decreases,  wiH  have  less  damping  effect. 

Therefore,  terms  for  large  values  of  x  cannot  be  neglected.  The  detectors  farther 
fran  the  center  will  weigh  more  heavily,  while  the  center  detector  wei^ting  co¬ 
efficients  or  detector  sensitivities  will  be  smaller.  On  the  other  hand,  if 
m 

—  is  small,  the  damping  effect  will  dominate,  and  a  few  detectors  near  the  center 
would  be  sufficient. 

From  Equation  (ll),  it  is  seen  that  the  inverse  transform  of  the 

band-pass  filter  is  a  "carrier"  of  frequency  2™^,  amplitude-modulated  by  a 

"signel"  of  mian.  The  majcimum  detector  width  usable  is  determined  by  the  center 

frequency,  and  is  equal  to  —  .  Larger  width  detectors  would  not  fit  b’oween  all 

™c 

the  zeros  of  the  filter  function.  Am,  the  bandwidth,  should  detennine  the  trunca¬ 
tion  point,  since  for  smaller  bandwidths,  the  cutoff  rate  becomes  increasingly  im¬ 
portant. 

III.  NUMERICAL  SOLUTION 

Computing  F*  (jm)  bj  hand  would  have  been  very  tedious.  Therefore,  an  IBM 
709^  was  used  o  make  a  numerical  computation.  Since  it  seemed  desirable  to  evalu¬ 
ate  the  effects  of  truncation  and  detector  width  separately,  a  doubly  parametic 
series  of  curves  was  chosen.  One  parameter  is  the  normalized  detector  width,  Wm^. 

The  second  parameter  is  the  normalized  truncation  poin1^ 

A.  LOW  PASS 

Assuming  a  low-pass  filter  with  a  cutoff  frequency,  m  ,  the  zeros  of 
f(x)  appear  at  x  =  ,  where  k  is  an  integer  ranging  from  1  to  od  .  Normalized  de¬ 

tector  widths  of  Wm^^qual  to  O.O5,  0.1,  0.25,  and  O.5  were  chosen  (see  Figure  2), 
and  also  normalized  truncation  points  equal  to  1.0,  2.0,  and  k.O.  Ilie  results 

ere  shown  in  Figures  4,5>  and  5. 

1.  Effect  of  Detestor  Width 

For  a  given  truncation  point,  for  example,  x^  »  1  (see  Figure 
'*),  the  approximation  is  recognizable,  but  not  good  when  oily  four  detectors  (2 
each  side,  Wm  ■  0.5)  used.  As  the  number  of  detectors  is  increased,  the  ap- 

w 

proxxmation  is  improved,  particularly  in  the  stop-band,  because  the  harmonic  maxima- 
rainima  are  moved  to  larger  veulues  of  m  where  they  are  attenuated.  However,  increasing 
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the  number  of  detectors  effects  a  diminishing  gain  in  pass-band  approximation  ac- 
cur-aey.  Again,  referring  to  Figiire  4,  the  improvement  is  considerable  when  the 
number  is  increased  from  4  to  8,  but  barely  discernible  when  increased  from  20  to 
40,  and  the  small  gain  hardly  justifies  the  numerical  increase  in  det  ctors.  Tnis 
effect  is  independent  of  truncation  point  as  seen  from  Figures  5  and  6. 

?.  Effect  of  Truncation  Point 

For  a  fixed  detector  width  (e.g.,  Wm^  =0.5),  an  increase  in  the 
truncation  point  such  as  shown  in  Figures  4,  5,  and  6,  causes  the  slope  at  the  cut¬ 
off  point  to  become  sharper  aind  sharper,  and  the  number  of  ripples  in  the  bass -band 
to  increase.  For  a  truncation  of  equal  to  only  1,  the  bass-band  response  is 

relatively  smooth. 

5.  Practical.  Application  of  Results 

In  selecting  the  proper  filter  characteristics  for  a  specific 
application,  both  the  pass-band  and  the  stop-band  requirements  must  be  specified. 

The  detector  size  wilT  result  from  the  permissible  behavior  in  both  the  pass-band 
and  the  stop-band.  If  the  first  harmonic  maxima -minima  must  be  a  given  number  of 
octaves  from  m^,  the  detector  size  is  chosen  accordingly.  Pass-band  ripple  magni¬ 
tude  will  also  establish  detector  size.  If  a  given  cutoff  rate  is  desired,  the 
truncation  point  is  established. 

B.  BAND  PASS 

Since  the  band-pass  filter  equation  is  derived  from  the  difference  of 

two  low-pass  equations,  the  band-pass  curves  can  be  very  nearly  predicted  from  the 

difference  of  two  scaled  low-pass  filter  cujr'/’es.  Four  normalized  detector  widths 

(0.25,  0.125,  0.05,  and  0.025)  were  chosen  for  the  numerical  band-pass  filter 

characteristic.  Two  bandpasses  -  a  narrow  band  =  0.2)  and  a  wide  beind  (—  =  l.O) 

c  “c 

Were  investigated.  Truncation  points  of  1,  2,  and  4  were  used  for  the  wide  band, 
while  values  of  5>  10,  and  20  were  used  for  the  narrow  band. 

1.  Effect  of  Detector  Width 

^rcxn  Figures  J  throu^  12  it  can  be  seen  for  both  the  wide  and 
the  narrow  bands  that  as  the  detector  width  decreases,  the  harmonic  maxima  and  minima 
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spread  and  reduce  in  ma^iitude,  in  a  fashion  similar  to  that  for  the  low-pass  cur\'es 

The  stop-hand  in  the  region  throu^  8  ^which  is  more  than  2  octaves  above  the 

cutoff  frequency)  for  the  smaller  detectors  wheie  is  equal  to  0.05  and  0.025, 

is  effectively  zero.  Above  8  however,  there  are  periodic  maxima  and  minima 

■vdiose  magnitudes  are  less  than  0.2  and  which  become  smaller  with  increasing  m. 

When  the  lower  resolution  detectors  are  used  (i.e.,  Wm  =  0.25  and  0.125  unit  length 

c 

the  approximation  is  very  sensitive  to  detector  width.  For  the  smaller  detectors 
■vAiere  Wm^  =0.05  and  0.025,  the  approximation  improvement  -  particularly  in  the 
pass  band  -  becomes  negligible.  Figures  J  throu^  12  all  illustrate  the  diminis]'.- 
ing  effect  of  decreasing  detector  width  for  fixed  truncation. 

2.  Effect  of  Truncetion  Point 

Here  also  the  results  are  what  mi^t  be  expected  from  extending 
the  low-pass  case.  Increasing  the  truncation  point  sharpens  the  edges  of  the  pass- 
band  characteristics  aind  also  increases  the  number  of  ripples  in  the  pass  band. 

For  the  maximum  detector  width  of  Wm^  =0.25  (see  the  wide-band  filter.  Figures  1, 

8,  and  9),  the  pass  band  exhibits  an  average  decreasing  magnitude  for  the  larger 
truncations,  though  not  for  the  minimum  truncation  of  =  !•  This  "slump"  in 

the  pass  band  is  corrected  by  using  smaller  detector  widths.  Truncation  and  de¬ 
tector  width  therefore  cannot  be  considered  as  having  entirely  independent  effects. 

5.  Practical  Application  of  Results 

Narrow-band  spatial  filters  would  be  very  useful  in  making  spatial 
frequency  measurements.  An  instrument  could  be  constructed  using  variable  weights 
for  each  detector  channel,  and  could  thereby  be  made  to  measure  a  set  of  frequency 
bands  in  turn.  1+  is  encouraging  that  even  for  minimum  truncation,  -  5  (see 

Figure  lO),  a  good  narrow  band  pass  can  be  obtained,  even  though  it  does  not  well 
approximate  the  square-sided  ideal.  For  larger  truncations  and  more  detectors,  the 
approximation  is  excellent. 

IV.  CONCLUSION 

The  synthesis  of  two  basic  one- dimens icaial  spatied  filters,  low-pass  and  band¬ 
pass,  using  an  array  of  detectors  has  been  demonstrated.  The  synthesis  is  truly 
spatial  in  .  at  there  is  no  functional  dependence  'n  time,  and  ♦'herefore  the  filter 
characteristic  is  independent  of  fluctuations  in  rate  of  scanning  mecheuiisms. 
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The  output  signal  of  a  spatial  filter  is  proportional  to  the  spectral  power  of  the 
observed  scene  within  the  pass  band  of  the  filter. 

Although  only  the  two  basic  filters  were  discussed,  finite  size  detectors 
make  it  possible  for  almost  any  finite  bandwidth  filter  characteristic  to  be  ap¬ 
proximated  and  synthesized  by  this  method.  A  specific  application  would  involve 
design  of  a  filter  to  best  discriminate  against  its  particular  background  distri¬ 
bution.  Any  such  general  filter  approximation  would  be  dependent  on  the  individual 
detector  field-of-view  and  also  on  the  field-of-view  f  the  array.  A  practical 
utilization  of  such  a  filter  would  amplify  each  output  separately  (because  of 
signal-to-noise  arguments),  weight  the  outputs,  and  sum. 

If,  by  scanning  a  detector,  the  spatial  function  is  described  as  a  parametric 
function  of  time,  then  the  filtering  can  be  theoretically  accomplished  with  a  set 
of  parallel  delay  lines  and  associated  gains  whose  outputs  are  summed.  This  method 
is  ccanpletely  equivalent  to  using  the  detector-array  spatial-filter  synthesis,  and 
results  from  the  same  principle. 
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FILTERS  WITH  CIRCUIAR  SYMMETRY 

I.  TWO- DIMENS  lONMi  CIRCUIAR  SYRg^ETRY  LOW -PASS  FILTER 

The  synthesis  of  a  two-dimensional  low-pass  filter  with  a  cutoff  frequency, 
me-,,  on  the  x  axis  and  a  cutoff  frequency,  Y  axis,  ■was  discussed  in 

Section  III.  The  cutoff  frequency  for  any  other  axis  is  the  intersection  of  tliat 
axis  and  the  cube. 

For  a  circularly  symmetric  low-pass  characteristic,  that  ■’s,  one  independent 
of  the  direction  of  scan,  the  response  function  is  a  cylinder  in  the  frequency 
domain.  The  cylinder  has  a  height,  1,  for  p<  where  p  is  the  radial  frequency 
in  any  direction  and  p^  is  the  cutoff  frequency.  See  Figure  8. 

The  inverse  transform  of  this  function  in  polar  form  is  f  (r,  0)  where 


,0 


f  (r,  9)  = 


0  dp 


O  O 


rP, 


'JT  p  (2tT  pr)  dp 


(B-1) 


-  f  Jj.  (2"  Pe  >■) 


(r,  0)  is  plotted  in  Figure  8. 
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Using  annular  detectors  of  -width  W  and  weight  ,  where  the  wei^^t  is  taken 


fr(Sn  the  value  of  f  (r,  9)  at  r  i  o  j 


h  '  Pc 


(B-2) 


The  inverse  transform  of  each  annular  detector  is 


F.  (p,  0)  = 


W 

(i-l)W 
2  P, 


^  P 

(n  p^  iW)  271  r  Jq  (-277  p  r)  dr 


J.  P.  iW)J 


p  “1 


(27t  p  iW)  -  (  jj27t  p  (i-l)w) 


(B-3) 


and 


F*  (p,  0) 


n  ^ 

=  ^  («  P^  iW)  (27T  p  iW)  -  J^(27r  p  (i-i)vj)>(B-4) 

i=l  L 


II.  TWO-DIMENSIQML  CIRCULAR  SYMMETRY  BAKD-FASS  FILTER 

For  an  annular  band-pass  between  p^  -  (Ap/2)  and  p^  +  (Ap/2)  the  inverse 
transform,  f  (r,  9),  is 


f(r,  e)  = 


o 

c  2  I 


2"  (po  +  t)p 


c  2 


(pc  ’  f)  ^ 


(B-5) 


.th 


With  annular  detectors  of  width,  W,  the  i  detector  has  a  weight,  A^,  of 


A.  - 
1 


Ap 

21Pc  +  Y 
iW 


2«(Pc  +  t) 


P  ^ 

iW 


"■'(pc  ■  f-) 


(b-6) 
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and  therefore 


F*  (p,  0  )  -  V \ 

ti 


U  •  f]  U  fcM 

iW  "^1  2  /  ^  iW 


/  ^ 

^1 


(B-7) 


/  \  r  1 

X  jjj^[2n  p  iw]  -  J^[2rt  p  (i-1)  . 

III.  TWO-EIMENSIQNAt.  BAND-P/^G  WITH  RECTANGULAR  SYJWTRY 

The  t’-ro-dimensional  hand-pass  space  function  is  the  difference  of  two  low- 
pass  space  functions.  That  is 

fr  /  .  V  n 


f  (x,  y)  =  ^ 

n 


sin  2n  m 


■of  ^  ,  4  o  (  ^  A 

\  =1  2  ;  1=2/ 


{b-8) 


The  weight  for  a  detector  whose  center  is  located  at  y^^  is 


A.,  = 

IK 


/  A  m  \  W  f  ^  % 

^  sin  2n  j  (2i-l)  sin  2n  U  +  ~^j  (2k-l)  ~ 

F  I  (21-1)  W^/2  (2k-l)  Wg/S 


sin  2s  — j  (21-1)  ^ 


W-J  /  ^  \  W 

(2i-l)  ~  sin  2n  (2k-l)  -f 


(2i-l)  W^/2 


(2k-l)  W„/. 


(B-9) 


and  the  trarisform  of  the  array,  F  (m^,  m^ )  is 
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F*  (m^,  rtip) 


1  sin  n  in,  \1 .  sin  n  W_  — ,  c-ii, 

L  _ l-A  ^  2  V  V 

2  7t  m  w  rt  W  L-i 

^  i=l  k=l 


X  cos  (2i-l)  cos  (2k-l)  (B-10) 
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APPENDIX  C 

RELATION  BETWEEN  MAXIMUM  SLOPE  OF  SPATIAL 
FREQUENCY  FUNCTIONS  AND  MAXIMUM  SCANNING  DETECTOR  DIMENSIONS 


We  already  know  that  the  spatial  fi^iquency  response  functions  of  scanning 
array  detectors  do  not  exactly  reproduce  the  desired  functions.  A  useful  param¬ 
eter  for  'Characterizing  quality  of  those  spatial  frequency  functions  that  are 
designed  to  have  band-  or  low-pass  characteristics  is  the  slope  of  F  (u)),  that  is 
dFliJ  , 

d  UD 

When  F  (u))  represents  the  Fourier  transform  of  an  arbitrary  space-limited 
function,  f(x),  an  interesting  relationship,  known  as  Bernstein's  theorem,  gives 
the  maximum  slope  of  F  (o))  in  terms  of  the  width  of  f(x).  Here  we  apply  this 
theorem  to  one-  and  two-dimensional  rectangular  coordinat"  Fourier  transforms, 
and  to  the  polar  coordinate  Fourier-Bessel  transform  (Hankel  Transform).  Its 
derivation  is  given  by  Arsac,  who  follows  the  work  of  Boas. 

For  a  one-dimensional  Fourier  transform  the  relationship  is 


d  F  (ui) 
d  ut) 


f  L  B  (F) 


(C-1) 


where  L  represents  the  half -width  of  »  function  f  (x)  that  extends  from  -L  to 
+L,  and  B  (f)  represents  the  maximum  value  of  I  F  (u))|  • 

__ 

J.  Arsac,  Transformation  de  Fourier  et  Theorle  des  Distributions 
(Paris,  Dunod,  I96I),  pp.  23^-259* 

R.  P.  Boas,  Entire  Functions  (New  York,  Academic  Press,  195^),  PP.  210,  211. 
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In  two-dimensional  rectangular  coordinates 


5  F  ot)^) 

^  U)  ^  (D 
X  y 


< 


L  L 
X  y 


B  (F) 


(C-2) 


where  L  and  L  represent  maximum  dimensions  in  x  and  y  directions. 

X  y 

For  a  circular  function  of  bounded  dexsain,  that  is,  a  filter  function  con- 

2  2  2 
fined  to  the  circle  x  +  y  <  R*" 


S  F  (u,^) 
^  (1) 

r 


<  JT  R  B  (F) 


(c-5) 


Although  the  slopes  of  the.  F  (uj)  functions  that  result  from  the  scanning 
weighted  detector  arrays  will  be  less  than  the  maxima  given  here,  the  bounds  can 
be  useful  in  guiding  preliminary  design. 
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COME  ARISON  OF  TRAI^PARENCIES  WITFI  ELECTRICALLY  WEIGHTED  DETECTOR  ARRAYS 

Since  most  of  this  report  uses  the  ideal  low-pass  spatial  frequency  filter 
as  an  example  for  specific  discussion,  we  shall  use  it  here  also.  Thus,  con¬ 
sider  the  function  f(x)  =  (l/it)(sin  x/x);  that  is,  the  inverse  Fourier  transform 
of  an  ideal  low-pass  filter,  shown  in  Figure  A -IB.  As  described  in  Appendix  A 
and  shown  in  Figure  A-2,  an  array  of  electrically  weighted  detectors  approximates 
f(x)  =  (l/jT)(sin  x/x)  by  a  number  of  rectangles,  each  representing  one  detector. 
The  height  of  a  rectangle  corresponds  to  the  electrical  weighting  of  its  detector 

arKi  can  be  positive  or  negative,  the  negative  weighting  realizable  by  electrical 
inversion. 

Instead  of  an  array  of  detectors  forming  f(x)  -  (l/n)(sin  x/x),  why  not 
cover  a  single  detector  with  a  screen  of  prescribed,  variable  transmission, 

T(x)  =  (sin  x/x;?  The  obvious  objection  is,  of  course,  the  fact  that  T(x)  can 
have  only  positive  transmission;  that  is,  0  <  T(x)  <  1.  As  we  show  here,  there 
are  ways  to  get  around  this  difficulty,  but  fJrst  let  us  point  out  the  limitation 
that  a  positive  function  imposes  on  its  Fourier  transform 

This  limitation  is  clearly  seen  by  examining  the  Fourier  transform.  Assume 
that  f(x)  has  been  normalized  so  that 

J+OB 

f(x)  dx  =  1 
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Tlien  frcsn 


F(u)) 


we  see  at  once  that  F(0)  =  1.  Next, 


(D-2) 


F(a;)  < 


'  +  00  C 

<  1 

sj 

f(x) 

dx  < 

f (x)  j  €  1  dx  < 

f(x)  dx  =  1  (D-5) 


using  the  fact  that  f{x)  is  real  and  positive.  Thus,  |f(u,)|<  F(0);  this  means 
that,  in  a  strict  sense,  F(ui)  cannot  form  a  high-pass  filter  function.  The 
qualification  is  needed  hecause  F(ui)  can  have  a  greater  magnitude  at  a  higher 
frequency  that  at  some  lower  one.  A  noteworthy  case  in  "vrhich  F(uj)  acts  as  a 
hand-pass  filter,  though  it  is  the  Fourier  transform  of  a  positive  function, 
arises  when 


f(x) 


1  -*•  cos  (ju 


-L  <  X  <  L 
X  L 


(D-  ) 


Its  transform 


’(o))  =  2  L 


sin  ^'L 


+  - 


sin  (a)-uiQ)L 


sin 

{x'+!JUq)L 


(D-5) 


amounts  to  a  pair  of  (sin  x/x)  functions  centered  at  ^  and  -aj  ,  and  a  larger 

0  0 

(sin  x/x)  centered  at  uj  =  0  with  twice  the  amplitude  of  the  pair  at  and  -uj^. 
Thus,  F((ji))  can  he  considered  as  a  hand -pass  function;  hut  the  large  amplitude 
pass -hand  at  (y  =  0  distinguishes  F(uj)  from  what  an  electrical  engineer  might 
conventionally  consider  as  a  hand-pass  filter. 


While  the  foregoing  example  does  indeed  show  how  a  positive  function  can, 
in  a  restricted  sense,  act  as  a  band-pass  filter,  it  is  in  fact  an  exceptional 
case.  Consequently,  we  return  to  the  discussion  of  how  to  realize  negative 
transmissi on. 
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*  -SHf 

Recent  work  ’  has  shown  some  practical  schemes  that  provide  negative 
transmission  in  image  processing  experiments.  But  this  work  is  concerned  with 
processing  in  which  the  negative  transmittance  effect  precedes  photographic 
recording  of  the  processed  image.  Our  need  for  negative  transmission  is  easier 
to  satisfy'  because  the  electrical  outputs  of  separate  detectors  can  be  subtracted 
Just  as  easily  as  added,  whereas  in  the  work  referenced  above,  subtraction  uses 
either  the  Herschel  effect  in  which  long  wavelength  light  erases  an  image  that 
has  previously  been  formed  with  shorter  wavelength  light  ;  or,  a  fluorescent 
screen  which  brightens  under  ultraviolet  radiation  and  darkens  under  infrared. 


With  two  electrical  detectors,  negative  transmission  can  be  achieved  in  the 
following  way.  Fii.®  t  (x,y),  the  desired  transmission  function,  form  the  two 
positive  functions  t  (x,y)  and  t_  (x,y),  where  t  (x,y)  =  t.  (x,y)  -  t  (x,y), 

‘  T  “ 

defined  as 


(x,y) 


1  t  (x,y)  where 
1  0 


t  (x,y)>  0 
elsewhere 


(d-6) 


(x,y) 


■rhere 


t  (x,y)<  0 
elsewhere 


(D-T) 


Next,  construct  transparencies  representing  these  functions.  Then  simply  split 
the  image  illumination  into  two  identical  beams;  send  one  beam  through  t_^  (x,y) 
onto  a  detector  whose  output  has  the  electrical  weight  of  plus  one;  send  the 
other  bea.m  throiogh  t_  (x,y)  onto  a  detector  whose  output  has  the  electrical 
weight  of  minus  one;  and,  finally,  s^im  the  two  detector  outputs.  This  process 
can  be  described  in  terms  of  the  Fourier  transforms  of  the  image  and  the  trans¬ 
parencies  as  follows. 


* 

1.  A.  Trabka  and  P.  G.  Roetling,  J.  Opt.  Soc.  Am. .  54,  1242  (1964). 
D.  H.  Kelly,  J.  Opt.  Soc.  Am..  1095  (I961). 
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Let  I  (iu  ,  0)  )  renresent  the  Fourier  transform  of  the  image  intensity.  Let 

T  (u)  ,  UJ  )  represent  the  Fourier  transform  of  the  transmission  function,  where 
X  'y 

T  -  T_  (u)^,  a)y)  (D-8) 

corresponding  to  t  (x,  y)  =  t  (x,  y)  -  t_  (x,  y).  With  0  (u)  ,  cu  )  the  processed 

+  ~  X  y 

transform,  wo  find 

0  (u)^,  uj^)  =  T  I  (ui^,  u)y)  =  i  -  T_  I 

This  equation  summarizes  the  negative  transmission  process. 

Having  described  one  way  of  realizing  negative  transmission,  we  now  discuss 
some  practical  differences  between  such  schemes  and  the  electrically  weighted 
detector  array.  First  of  all,  the  two -transparency  filter  reduces  detector 
illumination,  whereas  electrical  weighting  allows  each  detector  to  receive  full 
illumination.  In  those  applications  in  which  detector  and  preamplifier  (ampli¬ 
fication  prior  to  summation  of  detector  signals)  noise  is  significant,  the  reduced 
detector  signal  of  the  two -transparency  filter  presents  a  disadvantage  in  com¬ 
parison  with  the  electrically  weighted  detecto:  array.  (We  assume  equal  detector 
area  in  the  two  schemes,  but  factors  which  have  been  overlooked  here  may  make  this 
an  unreasonable  basis  for  comparison. ) 

Electrically  weighted  detector  arrays  offer  the  advantage  of  flexibility  in 
that  their  'weighting  can  be  varied  to  give  different  spatial  frequency  filter 
functions.  It  is  conceivable  that  such  filter  function  control  could  help  to 
correct  optical  non-isoplanaticism.  That  is,  if  a  filter  had  been  designed  to, 
in  some  prescribed  manner,  operate  on  an  optical  point  spread  function,  and  the 
spread  function  varied  frem  place  to  place  across  the  image,  then  an  electrically 
controllable  filter  function  could  compensate  such  variations. 

Another  advantage  of  electrically  weighted  detector  arrays  might  arise  if 
detectors  had  non-uniform  gain.  Thus,  electrical  weighting  could  compensate 
gain  variations  from  detector  to  detector  while  simultaneously  providing  tae 
weighting  demanded  by  the  transmission  function. 


(D-9) 
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Finally,  electrically  'weighted  detector  arrays  may  some  day  serve  in  image 
processing  systems  that  require  filters  with  transmission  functions  that  can  he 
altered  in  response  to  changes  in  Input  images.  The  disadvantage  in  the  inflexi¬ 
bility  of  the  two “transparency  fitter  is  evident  here,  as  well  as  in  the  two 
preceding  applications . 
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